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It is argued that the electron stripes as found in correlated oxides have to do with an unrecognized 
form of order. The manifestation of this order is the robust property that the charge stripes are 
at the same time anti-phase boundaries in the spin system. We demonstrate that the quantity 
which is ordering is sublattice parity, referring to the geometric property of a bipartite lattice that 
it can be subdivided in two sublattices in two different ways. Re-interpreting standard results of 
one dimensional physics, we demonstrate that the same order is responsible for the phenomenon of 
spin-charge separation in strongly interacting one dimensional electron systems. In fact, the stripe 
phases can be seen from this perspective as the precise generalization of the Luttinger liquid to 
higher dimensions. Most of this paper is devoted to a detailed exposition of the mean-field theory 
of sublattice parity order in 2-1-1 dimensions. Although the quantum-dynamics of the spin- and 
charge degrees of freedom is fully taken into account, a perfect sublattice parity order is imposed. 
Due to novel order-out-of-disorder physics, the sublattice parity order gives rise to full stripe order 
at long wavelength. This adds further credibility to the notion that stripes find their origin in the 
microscopic quantum fluctuations and it suggests a novel viewpoint on the relationship between 
stripes and high Tc superconductivity. 
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I. INTRODUCTION. 

This rather long paper is devoted entirely to a possi- 
ble answer to the question what are stripes ? This might 
sound odd since the electron stripes as found in corre- 
lated oxides have grown into a popular subject. Of course 
there is a popular answer to this question: the 'rivers of 
charge' separated by (quasi) insulating domains (Zaanen, 
1999 (2)). However, this answer is too intuitive and lacks 
the precision needed in the context which matters most: 
are stripes central to the problem of high Tc supercon- 
ductivity, or do they represent a red herring, a real but 
irrelevant side effect? 

In the light of the history of the subject, nobody can 
afford to be convinced of anything related to high Tc su- 
perconductivity. Nevertheless, one might want to stick to 
the most general principles, namely those of symmetry. 
The high Tc riddle has to do with the highly anoma- 
lous macroscopic properties of the electron system of the 
cuprates. Hence, the question is on the long wavelength 
properties of the system and here one is helped by gen- 
eral symmetry-based considerations, of the kind called 
by Laughlin 'competing orders' (Laughlin, 1998). There 
appear to be two possibilities: either the fixed point is 
adiabatically connected to the ECS superconductor or 
the fixed point is a different one. In the first case, one 
has to explain why the fixed point is approached by a 
highly anomalous cross-over and this suggests the prox- 
imity of a different fixed point. In the second case, high 



Tc superconductivity is separated from a conventional 
superconductor by a non-adiabatic boundary and this 
implies that the symmetries of both states are different. 
The high Tc superconductor is surely a Meissner phase, 
and the difference in symmetry is somewhere else. It is 
apparently hard to detect by experiment and therefore it 
is called 'hidden order'. Also if it is a mere crossover be- 
havior, the nearby competing state is apparently equally 
hard to detect experimentally and one might want to 
consider it as a variation on the hidden order theme. 

Let us consider stripes from this perspective. Are they 
candidates for the hidden order? It is a popular thought 
that stripes are just charge order (or charge density wave 
order, or just an interesting Wigner crystal: it is all 
the same). The charge density of the electrons breaks 
translational invariance. Given the stability of the Mott- 
insulator it is not remarkable that this charge density 
wave is such that the hole poor regions become charge 
commensurate with the crystal lattice. Hence, these turn 
into magnetic domains and the subsequent breaking of 
spin-rotational invariance is then considered to be a par- 
asitic effect. 

Can this charge order be the hidden order? Despite 
some brave suggestions (Castellani, Grilli and di Cas- 
tro, 1995), this appears as highly unlikely. Translational 
symmetry breaking is easy to detect and it is not seen 
in experiment in the vicinity of optimal doping in the 
best superconductors. Of course, this could have been 
missed in the scattering experiments for technical rea- 
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sons. However, it is by now well established that nodal 
fermions exist also at lower dopings (Orgad et at, 2001) 
and this fact is very hard to reconcile with a significantly 
developed stripe charge order. 

An alternatively candidate could be the stripe mag- 
netism. In fact, the incommensurate magnetic fluctua- 
tions (often associated with the stripe magnetism) be- 
have in a manner which is reminiscent of the compet- 
ing order hypothesis. The case has been made that they 
demonstrate the 'quantum critical' scaling behavior asso- 
ciated with the proximity to a continuous quantum phase 
transition (Aeppli et al, 1997). However, to arrive at a 
more detailed interpretation one has to invoke a strong 
spin-charge separation (Chubukov and Sachdev, 1993; 
Sachdev, 2000; Zaanen, 1999(1); van Duin and Zaanen, 
2000) . The magnetism goes its own way (presumably de- 
scribed by a quantum non-linear sigma model) regardless 
of the charge dynamics. This hypothesis is directly vio- 
lated by experiment. It is clear that the gap seen in the 
spectrum of incommensurate spin fluctuations opens up 
right at the superconducting transition (Dai et al, 1999; 
Lee et al., 2000). 

Hence, it can be argued successfully that neither the 
charge order, nor the magnetic order associated with the 
stripes can be held responsible for the long wavelength 
anomalies of the high Tc phenomenon. The reason is that 
the empirical consequences of these conventional orders 
are too well understood and too easy to measure. 

The main aim of this paper is to illustrate the idea that 
the above might be an incomplete characterization of the 
symmetry structure of the stripe phase. Stripe order im- 
plies that yet another symmetry is spontaneously broken 
and this symmetry structure is of a most unconventional 
kind. On a heuristic level it is widely recognized that 
something unusual is going on and this is called 'topo- 
logical doping' or 'anti-phase boundarieness' (Zaanen and 
Gminarsson, 1989: Kivelson and Emery, 1996; Zaanen, 
Horbach and van Saarloos, 1996; Zaanen, 1998; Pryadko 
et al., 1999). All available experimental information sup- 
ports the notion that the charge stripes are at the same 
time anti-phase domain walls in the anti-ferromagnet 
(Zaanen, 2000(2)). This anti-phase boundarieness is ro- 
biist. It is not only there for static stripes (Tranquada et 
al., 1995, 1999; Emery, Kivelson and Tranquada, 1999); 
the 'dynamical stripes' are also defined through this anti- 
phase boundarieness. Most of the information on the 
latter comes from the spin-fluctuations as measured by 
inelastic neutron scattering. The interpretation of these 
in terms of stripes rests in flrst instance on the character- 
istic wave-vectors of these fluctuations and their depen- 
dence on doping (Aeppli et al, 1997; Mook and Dogan, 
1999; Mook et al. 2000; Dai et al. 1999; Yamada et al, 
1998; Lee et al, 2000). This assumes anti-phase bound- 
arieness which extents up to large energy scales. This 
anti-phase boundarieness is also a common denominator 
in many theoretical works addressing the microscopy of 



the stripe phase, ranging from the early mean-field work 
to the sophisticated recent exact diagonalization stud- 
ies (Zaanen and Gunnarsson, 1989; Zaanen, 1998; White 
and Scalapino, 1998; Morais-Smith et al. 1998; Pryadko 
et al, 1999; Fleck et al, 2000; Martin et al, 2000; Sto- 
jkovic et al., 2000; Tchernyshyov and Pryadko, 2000). 
Here we will largely ignore stripe microscopy and instead 
try to contribute to the understanding of the long wave- 
length dynamics from a phenomenological perspective. 
We want to suggest that this anti-phase boundarieness is 
a manifestation of spontaneous symmetry breaking. 

Stripes are not completely on their own in this re- 
gard: the same symmetry principle is behind the phe- 
nomenon of spin-charge separation in the Luttinger liq- 
uids of one dimensional physics. In fact, symmetry-wise 
stripes should be understood as the unique realization of 
spinful Luttinger liquids in 2-l-lD. 

The expert in one dimensional physics might find this 
puzzling: why should spin-charge separation have any- 
thing to do with spontaneous symmetry breaking? This 
will be discussed in some detail in section II. We start out 
with well-known, mathematical results in 1-l-lD physics 
(Ogata and Shiba, 1990) to find that these can be refor- 
mulated in the language of order parameter correlators 
after identifying the degree of freedom which is carry- 
ing the order. This degree of freedom is quite simple 
but counter-intuitive: for lattice problems it is sublattice 
parity (Zaanen and Nussinov, 2001). 

Hence, we claim that the Luttinger liquid is charac- 
terized by an order parameter. Sublattice parity is an 
Ising (Z2) degree of freedom and true long range order 
can therefore exist at zero temperature in 1-1- ID. In the 
order parameter formulation it becomes trivial to gener- 
alize spin-charge separation to higher dimensions and in 
2-l-lD this turns into the anti-phase boundarieness of the 
stripes. 

As will be discussed in section II, one might want to 
view this sublattice parity order theory as a parametriza- 
tion in terms of auxiliary degrees of freedom of a more 
hmdamental geometric structure. 'Geometric' is used 
here in the same sense as in the Einstein theory of grav- 
ity. The geometry of embedding space is different for 
an external observer and the internal observer, in the 
present case the experimentalist and the spin system, re- 
spectively. However, this geometry is very simple, at least 
as compared to that of fundamental space-time. The spin 
system is a quantum antiferromagnet and the only prop- 
erty of embedding space which matters is the bipartite- 
ness of the effective lattice seen by the spins. Full order in 
the sublattice parity language implies that the true lat- 
tice seen by the spin system is bipartite everywhere and 
thereby geometrically unfrustrated. In one dimension it 
is always possible to divide a lattice in two sublattices 
and we suspect that this offers a special protection to the 
spin-charge separation phenomenon. However, in higher 
dimensions bipartiteness is not automatic. In the order 
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parameter language, topological excitations of sublattice 
bipartiteness order can be identified in 2+lD which can- 
not exist in I+ID. These correspond in the geometrical 
language with curvature events, equivalent to essential 
frustrations in the spin system (Zaanen and Nussinov, 
2001). Although it demonstrates that in 2+lD sublat- 
tice parity order is not generic, the precise nature of the 
disorder theory is at present very poorly understood - it 
is a most unusual structure. We suspect that this struc- 
ture has something to do with the mysteries of high Tc 
superconductivity, and in the final section we will discuss 
some work in progress to illustrate the problem. 

The remainder of this paper (Sections III- VII) might 
be considered as a review on the part of the problem we 
understand fairly well. It summarizes a large amount of 
work carried out during the last 7 years by our group in 
Leiden (Zaanen, Horbach and van Saarloos; 1996; Zaa- 
nen, 1996; Eskes et ai, 1996; Zaanen and van Saarloos, 
1997; van Duin and Zaanen, 1998; Eskes et al., 1998; 
Zaanen, Osman and van Saarloos, 1998; Tworzydlo et 
ai, 1999; Zaanen, 2000) When we started this pursuit 
we were convinced that we were addressing rather unre- 
lated parts of the physics. Amusingly, as a lucky circum- 
stance we just looked at all the bits and pieces needed 
to arrive at a synthesis which we recognized only quite 
recently. It starts out with a single assumption, defin- 
ing an important limiting case of the general problem: 
in sections III- VII we present the mean-field theory of 
sublattice parity order in 2-l-lD. All we require is that 
sublattice parity order is perfectly obeyed in 2+lD. We 
even do not allow for local violations of this order, let 
alone the global violations as discussed in the previous 
paragraph. At first sight this amounts to a quite mild, 
partial constraint on the dynamics. The charges have 
still the freedom to delocalize and the spin system is a 
highly quantum-mechanical entity as well. The sublat- 
tice parity order just amounts to the requirement that 
the charges have to form connected lines, like Ising do- 
main walls. A rather lively 'toy universe' emerges and 
we refer to section III for an overview. Nevertheless, the 
sublattice parity order exerts an unexpected dominance. 
The mere presence of perfect sublattice parity order forces 
under all circumstances the charge and the spin to order 
as well. This is due to a rather counter-intuitive order- 
out-of-disorder mechanism: the more severe the micro- 
scopic quantum fluctuations, the more robust the long 
range order at large scales. Although this story does not 
solve any of the pressing problems in the high Tc context 
we do find it interesting. Albeit being the wrong limit, 
it corresponds with a reference point for the construc- 
tion of a more complete theory which might relate to 
high Tc superconductivity and we find the complexity of 
this mean-field theory fascinating. On a practical level, 
it adds a counter-intuitive meaning to the notion that 
stripe long range order can originate in the microscopic 
quantum fluctuations of spins and holes. 



II. SUBLATTICE PARITY ORDER. 

In this section we will not come up with anything tru- 
ely new. It is just a recollection of some well known facts 
of one dimensional physics. However, our consideration 
is focused on exposing the bare essence of what is meant 
by charge-spin separation. This will turn out to be a re- 
markably simple principle which can be trivially imposed 
in higher dimensions. A more extensive discussion will 
be pubhshed elsewhere (Kruis, 2001). 

Charge-spin separation refers to the general property 
of one dimensional electron systems that the electron is 
an unstable particle while at the same time the real prop- 
agating excitations of the system carry fractions of the 
charge of the electron and the spin separately (Ander- 
son, 1997). This property can be deduced in various 
ways. It was first found in the bosonization framework 
and later confirmed in the exact Bethe Ansatz solutions 
(Voit, 1994). It is clearly related to a topological struc- 
ture. Remarkably, this has to do with a kink topology, the 
topological structure associated with an Ising type field. 
This should not be considered as self-evident; the mani- 
fest symmetry of the problem is SU{2) x U{1) (the spin- 
and charge global symmetries), and why should there 
be a Z2 Ising topological structure at work? Exactly the 
same problem is encountered in the stripe context in 2-1-1 
dimensions. For every theoretical purpose, the spin sys- 
tem in the cuprates is SU{2) invariant. Why is it so that 
stripes are like Ising domain walls? 

In the one dimensional context one gets a first glimpse 
of the answer by considering a very simple and well 
known example. Consider a ID antiferromagnetic chain 
of Ising spins and remove a single spin somewhere as de- 
picted in Fig. (|l]). Now consider what happens when the 
spin vacancy or hole moves to the left. The spin neigh- 
boring the hole moves in the opposite direction and after 
a couple of hops one finds the hole surrounded by anti- 
parallel spins while two parallel spins reside at the origin. 
Another way of calling this fact is that the electron has 
split apart in a pure Sz = 1/2 excitation (the spin do- 
main wall at the origin, or 'spinon') and a pure charge 
excitation (the 'holon') because the spin of the original 
electron and that of the spin domain wall carried around 
by the hole (the anti-parallel spins surrounding the hole) 
compensate. One directly infers that this has nothing 
to do with the Ising symmetry of the spins. All what 
matters is that the spins have anti-ferromagnetic inter- 
actions while the motion of the vacancy is accompanied 
by the backward motion of the spin. The net effect of at- 
taching the domain wall to the hole is that the backward 
moving spin ends up in the same, predominantly antifer- 
romagnetic orientation relative to its neighbors after the 
hole has passed. Insisting that the spins surrounding the 
hole are parallel would imply that every move of the hole 
would shake-off a spinon and this would be a very costly 
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FIG. 1. Cartoon picture of the mechanism of spin-charge 
separation in one dimension. A hole is injected in an antifer- 
romagnetic S = 1/2 spin-chain (top). When the hole moves 
to the right the spin moves backward (middle). The result of 
this kinematical process is that a spin-domain wall carrying a 
5 = 1/2 quantum (spinon) is left at the origin while the hole 
is bound to a spin anti-domain wall, and this composite only 
carries charge (bottom). 



affair. 

The implication is that this simple effect also applies to 
Heisenberg spins. Imagining, for whatever reason, that 
these spins could be made to order in a perfect antiferro- 
magnet, we could have as well chosen to orient the spins 
along, say, the x direction instead of the z direction and 
the hole would still have been surrounded by anti-parallel 
spins. It is a less trivial matter to see that this effect is 
robust against quantum spin fluctuations. This has been 
proven rigorously to be the case, as we will discuss later, 
but it can already be inferred using a simple continu- 
ation argument for the single hole case. Although the 
microscopic spin fluctuations for the S = 1/2 case are 
severe, at long distances the classical Neel state is closely 
approached as signaled by the algebraic spin-correlations 
(Haldane, 1981). The hole domain wall is topological so 
that it exerts its influence at infinity, and at infinity the 
spins are closely approaching the classical limit. Con- 
sider the following Gedankenexperiment. Start out with 
the Ising spin chain and inject a single, completely de- 
localized Ising-holon which is, however, constrained to 
move in between two charge-potential barriers placed far 
apart. By measuring the spin correlator at two points 
outside and at opposite sides of these barriers one eas- 
ily infers that a holon has to be around because an up 
spin resides in the region to the right on the down spin 
sublattice of the region to the left and vice versa. Switch 
on subsequently the XY terms in the spin Hamiltonian. 
Upon reaching the Heisenberg point the spin correlations 



change from true long range order to algebraic order. It 
has to be that the kink is still around because it violates 
the algebraic correlations in the same way as it violates 
the true long range order of the Ising spins - a single Ising 
kink in one dimension suffices to cause true disorder. 

The nature of the Ising field supporting the kink re- 
mains to be clarified - from the previous discussion it is 
clear that this Ising field is unrelated to the spins them- 
selves. Although implicit to considerations of the kind 
discussed in the previous paragraph (Schulz, 1993), it 
seems that this field has not been explicitly identified be- 
fore: it is sublattice parity. Sublattice parity is defined 
as follows. Subdivide the lattice in A and B sublattices. 
Take an arbitrary reference point and start counting with 

either A - B - A - B ovB-A-B-A--- - and 

call the two possibilities 1 and —1, respectively. This is 
an Ising variable, except for the subtlety that the global 
degeneracy is a gauge degeneracy associated with the ar- 
bitrariness of the reference site: one could have as well 
started the counting from a neighboring site. The holon 
is a hole bound to a domain wall in the sublattice par- 
ity: the A sublattice changes into the B sublattice upon 
traversing the hole and vice versa. The simple kinemati- 
cal effect discussed in the previous paragraph translates 
into a geometrical principle governing the collective dy- 
namics. The only property of the embedding space which 
matters for a lattice quantum antiferromagnet is its bi- 
partiteness. The charge of the electron 'curves' the space 
as seen by the spin system, because it flips the parity of 
this bipartiteness. 

A much better job can be done, and this was already 
accomplished some time ago by Ogata and Shiba (Ogata 
and Shiba, 1990; Anderson, 1997). They observed that 
the Bethe-Ansatz ground state wave function of the Hub- 
bard chain has the remarkable property in the large U 
limit that it factorizes in a spin and a charge part. Con- 
sider a Hilbert space spanned by Ising spin configurations 
and holes. In the large U limit, as a first step one assigns 
the position of the holes and according to the Bethe- 
Ansatz analysis this configuration has an amplitude equal 
to the amplitude of an equivalent system of hard core 
bosons, regardless the configurations of the spins. Keep- 
ing the holes fixed, the spin-amplitudes follow from a 
pure Heisenberg spin system which lives on a chain which 
is geometrically altered. This is the 'squeezed' Heisen- 
berg chain, derived from the original Hubbard chain by 
removing the holes, and with the holes the sites where 
they reside, reinserting an antiferromagnetic exchange 
coupling J between the sites which were on both sides of 
the hole in the original chain, see Fig. ||. The true spin 
dynamics is the one of the pure spin system living on the 
squeezed chain, thereby explaining why spin and charge 
go their independent ways. As can be inferred from Fig. 
^, this squeezing operation is precisely equivalent to at- 
taching sublattice parity flips to the holes in the Hubbard 
lattice, which is just decoding the true space in which the 
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spin system lives in the fake space of experimentalists ob- 
serving the full Hubbard chain. Fundamentally, the prin- 
ciple underlying spin-charge separation seems to be best 
understood as a geometric principle - the space in which 
the spin dynamics lives is different from the Hubbard 
chain. However, the geometry involved is exceedingly 
simple (bipartiteness) and it is trivially parametrized in 
terms of an order parameter theory. This is an Ising the- 
ory. Every spin domain in between two holons is repre- 
sented by an Ising variable taking a value ±1 coding the 
value of the sublattice parity. Every hole charge is an an- 
tiferromagnetic exchange interaction between these Ising 
spins. In the large U limit there are apparently no fluctu- 
ations at zero temperature and a perfect Ising order is es- 
tablished: • • • (+l)-0-(-l)-0-(+l)-0-(-l)-0 . 

Spin-charge separation is just an Ising antiferromagnet! 

The examples discussed in the above both refer to 
rather specific situations (strong coupling Hubbard mod- 
els), and it is a-priori unclear if the Ising order dis- 
cussed in the previous paragraph is general. To get any- 
where, what is needed are order parameter correlators 
which can be explicitly calculated. These are necessar- 
ily of an unconventional kind: sublattice parity order 
can only be measured using the spin system because it 
parametrizes a geometric property of the space in which 
the spins live. The strategy is as follows: by inspecting 
the strong coupling limit a non-local correlator can be de- 
duced which removes the sublattice flips attached to the 
moving charges. This correlator thereby measures the 
true spin-correlations living in the squeezed embedding 
space-time, isolating its spin-only character. By inspect- 
ing this correlator one can indirectly infer the presence 
or absence of the sublattice parity order. 

Define the staggered magnetization as usual as Mx — 
{—l)'^Sx {x is the site index). Define also the charge 
operator taking the values 0, 1, 2 for an empty, singly, 
and doubly occupied site. Take an arbitrary point xq on 
the chain and define the following non-local (topological 
spin) operator (Zaanen and van Saarloos, 1997), 



T 



Mx 



T(xo,x) = e -^==0 M[x) 



(1) 



where the second line is the corresponding expression in 
the continuum limit. Now consider the correlator 



Otop{\x-x'lxa) = {■^\fx,^xfxo,x'\^) 



(2) 



the meaning of the 

'charge string operator' exp «7r^^^^^(l — n^^) is that 
it adds a minus sign every time that a hole is passed on 
the trajectory xq — x. One infers immediately that this 
charge string operator removes the flips in the sublattice 
parity attached to the holes. Instead of the antiferromag- 
netic sublattice parity order seen by the standard spin 
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FIG. 2. (a) According to the Bethe-Ansatz in the large 
U limit (Ogata and Shiba, 1990) one can first distribute 
the charges of the electrons and these charge configurations 
carry an amplitude equal to that of an equivalent system of 
hard-core bosons, regardless the nature of the spin configura- 
tions. The spin dynamics is that of a pure Heisenberg spin 
chain, where however the chain is a different one than the 
original Hubbard chain, (b) This 'squeezed chain' is obtained 
by removing the holes and the sites where the holes reside, 
substituting an antiferromagnetic exchange bond for these re- 
moved sites. This squeezing operation is precisely equivalent 
to attaching sublattice parity fiips to the holes in the original 
Hubbard chain. 



correlator • • • (+1) - - (-1) - - (+1) • • • the topo- 
logical correlator sees a ferromagnetic sublattice order 
• • • (-1-1) - 0- (-1-1) -0 - (-1-1) • • • because of the additional 
sign picked up every time a hole is passed. The topolog- 
ical correlator is easily evaluated in strong coupling and 
instead of the usual result [K^ and Kp are the spin and 
the charge stiffness, while e = 2kF — 7r/a = Tr/uh, where 
n/i is the hole density), 



Ospin(^\x X I) 



= B. 



cos(ea;) 



It is found that 



Otop{\x - x'\,Xo) = B„ 



1 - nh 
\x — x'\^'' 



(3) 



(4) 



Except for the 'dilution factor' 1 — rih Eq. (^ coincides 
with the correlator of a pure spin chain! The direct spin 
correlations as measured by Eq. (^ decay more rapidly 
because they are sensitive to the antiphase-boundarieness 
attached to the charges which are disorder events for the 
spins. Since the charges exhibit algebraic order this adds 
only a simple additional algebraic decay to the spin corre- 
lations ~ \x—x'\~^''. The charge-string operator removes 
the anti-phase boundaries from the spin sector and Otop 
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FIG. 3. In two dimensions, the charge string correlator cor- 
responds with a line integral and the correlations also depend 
on the length of the path L, see text. 

measures the physical spin correlations as they exist in 
the squeezed chain. 

By analyzing the standard lore of one dimensional 
physics (Voit, 1994) we have arrived at the conclusion 
that at least in strongly interacting lattice systems in 
1+lD charge-spin separation is controlled by a hidden 
order parameter. This order parameter structure is in 
turn breaking an Ising symmetry and should therefore 
be regarded as rather robust because it is protected by a 
mass-gap ~ it is the only true long range order which can 
exist in the one dimensional Luttinger liquid! Let us now 
proceed in a phenomenological fashion. States of matter 
can only be rigorously defined through their symmetry 
structure and we define the Luttinger liquid in arbitrary 
number of dimensions as states of matter which exhibit 
the same sublattice parity order as the Luttinger liquid 
in one dimension. This we find semantically more pre- 
cise than the widespread habit of attempting to define 
the Luttinger liquid in higher dimensions through the 
nature of the excitations (Anderson, 1997). The atten- 
tive reader should already have realized the inescapable 
conclusion: static stripes are the genuine generalizations 
of the Luttinger liquid to higher dimensions! 

Let us consider the generalization of the sublattice par- 
ity order to higher dimensions in more detail. Curiously, 
it appears to be not possible in all cases. Unlike the 
one dimensional scenario, not all higher dimensional lat- 
tices can be partitioned into two sublattices. It is far 
from clear how to construct stripe states on tripartite, 
etecetera, lattices and we suspect some interesting con- 
nections between the geometric frustration of the spin 
system and the topological interplay of charge and spin. 
It is likely not coincidental that the stripes have only 
been found up to now in crystals characterized by bipar- 
tite lattices. 

On square (cuprates, nickelates) or cubic (managan- 
ites) lattices, we can define a sublattice parity Ising vari- 
able just as in the one dimensional case. The charge 
is attached to domain walls in the sublattice parity and 



Ising domain walls are lines in 2D, sheets in 3D, etcetera: 
the stripes. As in the one dimensional case, by comparing 
the topological- and direct spin correlators one can learn 
if the sublattice parity order is established. Stronger, 
in space dimensions larger than one the topological cor- 
relator becomes more powerful. One can establish the 
sublattice parity order by inspecting only the topologi- 
cal correlator. Consider the generalization of Eq.(|^) to 
higher dimensions, 

Otopi\S-x'\,L) = (vI/|Al(f)e"^-^''''""^'^'Al(f')l*) 

(5) 

the charge-string operator corresponds with a line inte- 
gral and the correlations no longer depend only on the 
distance between the endpoints x — x' but also on the 
length of the path L of the path over which the inte- 
gral is evaluated (see Fig. H) . True long range sublattice 
parity is established if the following condition is satisfied, 

^ hm Otop{\x-x'lL)^G{\x~x'\) (6) 

\x—x'\^L — *oo 

where G is only a function of the distance between the 
end points. According to the present understanding of 
the stripe phenomenon, this condition should be satisfied 
in the static stripes of cuprates and nickelates. In fact, 
Zachar's recent analysis (Zachar, 2000) on the nature of 
the stripe disorder as driven by quenched disorder can be 
taken as a direct evidence that the condition Eq. (^ is 
satisfied in the static stripe phases of the cuprates. 

III. THE MEAN-FIELD THEORY OF 
SUBLATTICE PARITY ORDER. 

In the theory of order a most useful theoretical device 
is the limit where the order is perfect. This perfect order 
is barely ever realized. However, as long as the violations 
of the order are only local, the physics at long distances 
is qualitatively identical to that of the fully ordered case. 
For most of the remainder of this paper we focus on the 
consequences of perfect sublattice order in 2-1-1 dimen- 
sions. Hence, we impose that sublattice parity order is 
perfect which is equivalent to the statement that the in- 
ternal space after the Ogata-Shiba squeeze is a bipartite 
2D lattice, which is in turn equivalent to the statement 
that the charges are attached to Ising domain walls in 
the sense that they have to form connected d — 1 dimen- 
sional manifolds in the embedding space with d space 
dimensions. 

The shear length of this paper makes already the point: 
this mean-field theory of sublattice order is quite an in- 
teresting affair. The reason is that the presence of sublat- 
tice parity long range order leaves much room for other 
physics to happen, and the 2-|-l dimensional case is in 
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this regard more interesting than the Luttinger liquid of 
1+lD. As introduction to the remainder of this paper, let 
us introduce the several physics problems which emerge 
after imposing the sublattice parity order: 

(a) In principle, the rf — 1 dimensional stripe manifolds 
can have arbitrary shapes and the question arises what 
happens in such a system of interacting 'branes'. Specif- 
ically, in 2+1 dimensions stripes are lines on the time- 
slice and allowing for the fluctuations this represents a 
problem of interacting quantum strings, a string quan- 
tum fluid in 2+1 dimensions. This has been studied in 
great detail (Morais-Smith et al, 1998; Dimashko et ai, 
1999; Hasselmann et ai, 1999; Chernyshev et ai, 2000; 
Tchernyshyov and Pryadko, 2000), especially so in Lei- 
den (Zaanen, Horbach and van Saarloos, 1996; Zaanen et 
ai, 1996; Eskes et ai, 1996, 1998; Zaanen, Osman and 
van Saarloos, 1998; Zaanen, 2000). and we have acquired 
quite some insight in the nature of this problem. We were 
surprised several times. A first surprise is that the mere 
presence of an underlying lattice makes the problem of 
a single fluctuating stripe quite tractable. In fact, as 
compared to the continuum string theory of high energy 
physics these lattice strings are rather uninteresting ob- 
jects: they either pin to the lattice or they renormalize 
in free strings, as will be discussed in Section IV. In ad- 
dition, some general statements can be made about the 
system of interacting strings. Given the condition of com- 
plete connectedness, a case can be made that the system 
of strings in the presence of any interaction will always 
order (Zaanen, 2000; Mukhin, van Saarloos and Zaanen, 
2001). Given the results for the single string, the problem 
is obvious: a single free string already exhibits algebraic 
translational order, and it is obvious that the tenden- 
cies towards full order will be strong in a system of such 
strings. A particularly subtle case is the one where the 
strings only communicate via a non-crossing (hard-core) 
condition. This can be seen as as the decompactified 
(to 2+1 D) version of the gas of non-interacting spinless 
fermions in 1+lD and the argument will be reviewed in 
Section V, demonstrating that even this string gas even- 
tually orders, due to an order-out-of-disorder mechanism. 

(b) Although the spin system is at least globally un- 
frustrated in the presence of sublattice parity order (the 
'squeezed' lattice is bipartite) the quantum-spin physics 
in the 2+ ID case is still quite rich. The reason is that 
the stripes 'slice' the spin system in 1+1 D ladders, and 
the quantum-magnetism of the static stripe system can 
be discussed in terms of coupled ladders (Tworzydlo et 
ai, 1999; Sachdev, 2000), as will be reviewed in Section 

VI. A next problem is, what happens when the stripes 
are themselves strongly quantum fluctuating? In Section 

VII, we will present the results of a quantum Monte-Carlo 
simulation on a model which both incorporates the fluc- 
tuating stripes of Sections IV and V and the quantum 



spin dynamics of section VI (Osman, 2000). This has 
not been published before and it might be considered as 
the most sophisticated stripe model studied up to now. 
Besides illustrating vividly the physics discussed in the 
previous sections, it also adds a next piece of order-out- 
of-disorder physics: if the microscopic quantum-stripe 
fluctuations are sufficiently strong the spin system re- 
invents the classical Neel order, even if the spins of the 
fully static stripe system are quantum disordered! 

(c) Finally, it was implicitly assumed in the previous 
section that one full electron charge binds to every do- 
main wall unit cell, corresponding with the filled stripes 
of the nickelates. However, starting from the more gen- 
eral notion of the sublattice parity order there is no need 
to limit oneself to this special case. As a generalization, 
one might also want to attach some fraction of the elec- 
tron charge to the stripe unit-length. This might be half 
an electron corresponding with the half-filled stripes of 
the underdoped cuprates or even an irrational fraction so 
that stripes would be undoubtedly internally charge com- 
pressible metals - the metallic stripes of Kivelson, Emery 
et. al. (Kivelson and Emery, 1996; Emery, Kivelson 
and Zachar, 1997; Kivelson, Fradkin and Emery, 1998; 
Emery, Kivelson and Tranquada, 1999; Carlson et al.., 
2000; see also, e.g., Castro-Neto and Guinea, 1998; Za- 
anen, Osman and van Saarloos, 1998; Voita, Zhang and 
Sachdev, 2000; Fleck et ai, 2000; Bosch, van Saarloos 
and Zaanen, 2001). This adds yet another dimension to 
the physics and we leave a further discussion to these 
authors. We emphasize, however, that there is a-priori 
no conflict between this approach and what is discussed 
here. In fact, in the Emery-Kivelson school of thought, 
the local one dimensionality enters as an assumption. 
Sublattice parity order offers a rational for this assump- 
tion. 

IV. STRIPES AS STRINGS 

Let us take the Ogata-Shiba geometrical squeezing for 
granted, but now in 2+lD. The requirement that the 
squeezed lattice is an unfrustrated, bipartite lattice puts 
strong constraints on the way the stripes can fluctuate: 
only configurations are allowed where the 'holes' form 
fully connected trajectories, while every pair of holes 
is either nearest-neighbor ('horizontal bonds') or next- 
nearest-neighbor ('diagonal bonds'), see Fig. ^. 

This is a quite restrictive constraint. However, we no- 
tice that it is imposed by the requirement that sublattice 
parity is fully ordered. Longer excursions of holes away 
from the stripe would cause frustrations in the spin sys- 
tem on the squeezed lattice, thereby violating the order. 
Although in the cuprate reality these are likely to be 
quite important, our goal here is to derive the mean- field 
theory, and for these purposes it is necessary to neglect 
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FIG. 4. (a) A charged domain wall separating spin do- 
mains of opposite AFM order parameter, (b) If domain walls 
are 'broken up' spin-frustrations emerge and these events are 
excluded in the mean- field theory of sublattice parity order, 
(c) It is still possible to fluctuate the shape of the stripe by 
constraining its trajectory to only nearest- and next-nearest 
neighbor links, (d) A typical configuration is the kink and the 
movements of these kinks are an important source of kinetic 
energy, (d) An example of an undirected string configuration, 
(e) A typical directed string configuration. 



these motions. In addition, there are indications that the 
long wavelength physics of the single stripe is relatively 
insensitive for these 'microscopic details'. Hence, start- 
ing with a strong coupling model (with regard to binding 
of holes to stripes) we derive a fixed point physics with 
a finite basin of attraction. 

Despite the prescription that holes have to be nearest- 
or next-nearest-neighbors it should be immediately ob- 
vious that the stripe has still much room to quantum- 
fluctuate. What is the nature of the problem? 

Because of the connectedness requirement, at every in- 
stant of time the holes have to form 1+lD manifolds, and 
a single stripe is therefore a quantized string. As a for- 
tunate circumstance, the physics of a single string of the 
type following from the squeezing requirement has been 
already studied in a great detail. Eskes et ai, 1996, 1998, 
introduced precisely this kind of string for the purpose 
of a model study of the stripe fluctuations. Strings are 
extended entities which can exist in different collective 
states. The theory of strings, either in the high energy 
context or the membranes of statistical physics, is a rich 
subject which is far from completely understood. How- 
ever, it has long been recognized (Polyakov, 1987) that 
strings behave in one regard very differently from parti- 
cles. In the quantum theory of particles, the continuum 
limit can be reached by defining the theory on a lattice, 
taking subsequently the limit of the lattice constant go- 
ing to zero. This is in general untrue for strings and the 



richness of high energy string theory is in last instance 
associated with the true continuum limit. In the stripe 
context, the problem starts out with tight binding elec- 
trons moving on the crystal lattice. Hence, whatever else 
happens, the UV is lattice regularized and this turns out 
to be a most important condition determining the long 
wavelength dynamics of the fluctuating stripe. String 
theory on a lattice is easy. It coincides with the sta- 
tistical physics of crystal surfaces which appears to be 
completely understood (van Beijeren and Nolten, 1987). 

Eskes et at specialized on the particular lattice strings 
which are of direct relevance in the present context: 
neighboring 'holes' are connected by nearest- or next- 
nearest links on the lattice. Subsequently, Morais-Smith 
et ai, 1998, ( see also Hasselmann et ai, 1999; Dimashko 
et ai, 1999) studied a model with much weaker con- 
straints finding the same fixed points at long-wavelength 
suggesting that the constraints are microscopic details 
which are unimportant for the universal long wavelength 
physics. Eskes et al. managed to enumerate all possible 
phases of these strings. Although there is a wealth of 
ordered or partially ordered phases, corresponding with 
strings localized in space due to lattice-pinning, the only 
delocalizing string phase corresponds with the Gaussian 
fixed point, associated with algebraic long range order. 
Hence, on the single stripe level there is already a strong 
tendency towards localization and little is needed to find 
true long range order in a system of interacting stripes. 
To make matters worse, Eskes et al. found that at least 
these particular lattice strings are generically exhibiting a 
rather unusual type of symmetry breaking: these strings 
acquire spontaneously a direction in space, even if they 
are quantum delocalized. This is like a nematic order: 
although translational invariance is restored, rotational 
invariance is broken. If the string starts at —oo at the 
'left' of the 2D plane, it will always and up at the 'right' 
boundary, and never at the 'upper' or 'lower' boundary. 
This 'dircctedness' order originates in an order-out-of- 
disorder mechanism: directedness lowers the kinetic en- 
ergy at short distances. The statement can be made that 
it always happens because it works only better when the 
fluctuations become more severe. 

Let us present the mathematical definition of the Eskes 
strings, as well as a summary of the results insofar as 
they pertain to the remainder of this paper. Consider an 
configuration of points on the 2D square lattice, spanning 
up a ID trajectory where every pair of points is connected 
by either a nearest- (nn) or next-nearest-neighbor (nnn) 
link (Fig. The coordinate of the ^-th point is {r/f,rif) 
and the string configurations can be projected out from 
the total space of these points by. 



\r]{strings)) = ]JP(?7i+i - Vi) \v) , 



(7) 



where 
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(8) 



ensuring that the neighboring points are not farther apart 
than 1 or \/2 lattice constants. The potential energy of 
the string can be parametrized by, 



I 



1) 



(9) 



expressing that a nnn-link has an energy K relative to 
a nearest-neighbor one, and the lattice representation of 
curvature energy, 

2 

(10) 

expressing that e.g. two neighboring nn-links pointing in 
the same direction have a different energy than a nnn-link 
following a nn-link, or for instance two nn-links pointing 
in orthogonal directions. In principle one could also in- 
clude longer range link-link interactions but this will not 
change matters qualitatively at long wavelength. The 
string kinetic energy is, 



w2'' = rE^5t.(0^lt.(0 (e^*' 



(11) 



where tt is the canonical, periodic lattice momentum as- 
sociated with the position operator ff, 



(12) 



Acting once with T:f on a string configuration will cause 
a hop of point I over a lattice spacing in the x-direction, 
as long as the string constraint is not violated. 

This model is non-integrable and one can proceed 
in different fashions. In the path-integral formalism, a 
quantum particle corresponds with a worldline in a one- 
higher dimensional space, and likewise a quantum string 
becomes a worldsheet, a statistical physics membrane liv- 
ing in 2+1 dimensional embedding space. Lattice strings 
correspond with special membranes, namely those which 
also describe the statistical physics of crystal surfaces. 
The role of the lattice in the quantum problem is taken 
by the corrugation of the crystal in the crystal surface 
problem. 

It is easily seen that the general form of the action 
of the lattice string defined in the above is that of a 
restricted Solid- on- Solid (RSOS) surface problem. Here 
the surface is subdivided in columns with height rji and 
these column heights interact via terms hkc Eq.'s (|jl^) 
expressing that it costs for instance an energy K to have 
neighboring columns to differ in height by one unit, in- 
stead of having a flat configuration. It is also not hard to 
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FIG. 5. Illustration of the fact that a bend blocks the 
propagation of kinks along the string. Note that the ' holes' 
1 and 2 adjacent to the bend cannot move. 



find out that the lattice kinetic energy Eq. (|ll|) acquires 
a similar RSOS form after spreading it out along the time 
direction. 

A specialty of the lattice string is, however, that the t]^ 
and r]^ problems are described separately by their own 
RSOS surface and the interplay of the motions along the 
X and y directions gives rise to strong interactions be- 
tween both RSOS 'sectors' via local constraints. For in- 
stance, keeping both surfaces flat amounts to putting all 
particles I on the same lattice site. This problem was 
studied numerically, using quantum Monte-Carlo, and it 
, was discovered that in the parameter regime of interest 
always directedness symmetry breaking occurs. A par- 
ticularly interesting physical picture emerges in the lan- 
guage of coupled RSOS surfaces. In order to optimize 
the freedom to fluctuate, the best the system can do is 
to order one of the surfaces. In doing so, the constraints 
coming from the surface-surface coupling disappear com- 
pletely and the other surface can fluctuate freely. The en- 
tropy gained by this freely fluctuating surface out- weights 
the entropy associated with having both surfaces disor- 
dered. Take the x surface to be the ordered one. The 
order is such that this surface always steps upward, cor- 
responding with the string being directed along the x di- 
rection. Along the y direction the string can now freely 
quantum meander. 

The (directed) string can be viewed as a generalized 
quantum sine-Gordon problem, and it is most useful to 
consider its physics in terms of its soliton- or kink degrees 
of freedom. These are the events as shown in Fig. where 
the string steps sideways. It is easily checked that under 
the influence of the kinetic energy Eq. (^l|) these kinks 
propagate like free particles and this a dominant source 
of kinetic energy. To gain some intuition in the directed- 
ness symmetry breaking it is instructive to consider what 
happens when such a propagating kink approaches an 
'overhang' in the string, violating the directedness (Fig. 
H). It is easily seen that such an overhang acts like a 
hard wall for the soliton. This costs kinetic energy and 
this can be taken as an alternative physical picture for 
the mechanism driving the directedness. 
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At the same time, this directedness amounts to a great 
simphfication. The problem is reduced to a single RSOS 
problem and there is a great body of knowledge on RSOS- 
type models. It can be demonstrated that there are in 
total 10 distinguishable phases, see table (Q). Pending pa- 
rameters the string can be in various phases dominated 
by the potential energy where the stripes are localized in 
space. E.g., the string can be, on average, a straight line, 
which is pinned by the lattice, oriented along the hori- 
zontal (phase 11), or along a diagonal (phase I) direction 
in the lattice. However, also partially ordered phases are 
possible ('Haldane', 'Slanted' phases) and, last but not 
least, there is only one delocalized phase which is Gaus- 
sian as stated earlier. 

TABLE I. A schematic representation of the 10 different 
phases of the directed lattice string of Eskes et al. . Both char- 
acteristic configuration of the strings and that of the equiva- 
lent S = 1 chain are indicated. 



Phase String Spin 1 

I .• ++++++++ 

II 00000000 

III + + + + 

IV +0-+0+0-+ 

V -+0-0 0+- 

VI +-++-+ — 

VII . • • ' 0+0++0+0 

VIII . . • • " 0-FO+O+O+O 

IX 0+0+0-0+0 

X 0+0-0+0-0 



To get more insight in this phase diagram, it is in- 
structive to consider yet another representation of the 
problem: the directed string corresponds with a = 1 
Heisenberg spin chain with added Ising and single site 
anisotropics. This is easily seen in terms of a repre- 
sentation where the links are the dynamical degrees of 
freedom. Single out a particular 'guiding point' rjQ on 
the directed string and it is immediately clear that the 
string dynamics can be completely parametrized in terms 
of its center of mass 770 and the relative coordinates cor- 
responding with the set of link variables taking the val- 
ues 1,0,-1 corresponding with (l,f), (1,0) and (f,-l) 
bonds, respectively, for a string directed in the x direc- 
tion. For an infinitely long string the center of mass 
coordinate becomes non-dynamical, and the problem is 



completely parametrized in terms of the possible states 
on the links. These can be as well viewed as the three 
A/5 = 1,0,-1 states of a 5 = f quantum spin. For in- 
stance, the string kinetic energy is equivalent to the XY 
term in the S* = 1 spin representation, ~ S'j^ S^-^ + h.c. 

because S+ = V2(|I)(0| + |0)(-1|) in the basis of eigen- 
states of microscopic spin. Hence, acting once with this 
term changes two horizontal links into the two diagonal 
links corresponding with the sideward motion of the hole 
in the middle. 

The famous Haldane phase of the Heisenberg S = 1 
spin chain has a particular simple interpretation in the 
string language where it corresponds with a form of par- 
tial order (den Nijs and Rommelse, f989). In this phase 
(V in table I) , kinks have proliferated in the ground state 
and in this regard the state is quantum disordered. How- 
ever, there is still a form of hidden order in the sense 
that at average every kink which is moving the string 
upward is followed by a kink which is moving the string 
downward. Hence, the string as a whole is still local- 
ized in space although it is now localized in the mid- 
dle of two neighboring rows of the lattice (like a 'bond- 
ordered' stripe). This type of order is hidden from the 
spin-correlators and to make it visible in the spin chain 
one needs a non-local correlator. Eskes et al. discov- 
ered also a second type of partially ordered string: the 
'slanted string' (VII). This is like the Haldane phase ex- 
cept that the kink 'flavors' are now ferromagnetically or- 
dered such that the string orders along an arbitrary di- 
rection in the lattice. It was recently suggested that such 
a phenomenon might be relevant in the cuprate context 
(Bosch, van Saarloos and Zaanen, 2001). 

Most importantly, it is well established that S = 1 
quantum spin chains have only a single massless fixed 
point whose basin of attraction includes the XY point 
(phase IV) , where the string only has kinetic energy. This 
is a Gaussian fixed point and this is the only phase where 
the string is delocalized in space. This is an exceed- 
ingly simple fixed point: at large distances, the motions 
of the string can be completely parametrized in terms 
of the non-interacting transversal phonon-modes of the 
string. The position of points on the string can be writ- 
ten as r]{l) = ?7o(0 + Sr/{1), where 770 corresponds with 
the position of a flat string while Sr] corresponds with 
the transversal displacement. Following the standard 
lore of Gaussian theory it follows that the displacement 
correlator diverges logarithmically {{Sr]{l) — ^77(0))^) ~ 
ln{l) such that the string density correlator decays alge- 
braically (p(/)p(0)) ~ l/l^. 

The conclusion is, that a single string is either ordered 
or algebraically ordered. If one insists that stripes are 
fully connected entities, this might well be a very general 
conclusion. A major limitation of the above work is that 
it assumes the stripe to be internally incompressible and 
it is a-priori unclear what happens when the stripes are 
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internally metallic. A specific string model has been con- 
structed for a metallic stripe by ourselves (Zaanen, Os- 
man and van Saarloos, 1998). Although this also shows 
algebraic order at best, specific assumptions are made in 
its construction which render it to be less general and it 
might not be representative for the general case. 

However, on general terms it is not easy to see how 
one can avoid the algebraic string order. A sufficient 
condition for algebraic order is the non-vanishing of the 
line tension. How to get rid of a term in the action 
which is proportional to the world sheet area? There 
is actually one possibility which deserves a further explo- 
ration (Mukhin, van Saarloos and Zaanen, 2001): let the 
stripes be in equilibrium with a quantum gas of charges, 
while these charges can freely enter and leave the stripes. 
Under these conditions stripes are described by a qual- 
itatively different type of strings: the so-called extrinsic 
curvature strings. These are well known in the statis- 
tical physics context. Consider for instance biological 
membranes. These membranes are immersed in a 'gas' 
of constituents (lipids) which is in equilibrium with the 
membranes. If one pulls the membrane it simply absorbs 
lipids without paying a free energy penalty and therefore 
the tension (proportional to the membrane area) van- 
ishes. Instead, the next order invariant takes over, cor- 
responding with extrinsic curvature. After linearization, 
the action S ~ (i^^??)^, instead of the finite tension case 
S ~ {dfj,ri)^, and from power counting one infers directly 
that such an extrinsic curvature membrane fluctuates in 
the same way as a line with tension (like a worldline). 
The analogy with the stripe immersed in a bath of free 
charges should be immediately clear. There is however a 
caveat. The bare particles are non-relativistic and their 
action contains a mass term ^ {dr'rff' . Since the stripes 
are supposed to be made out of these particles it is hard 
to see how one can get rid of this tension in the time 
direction, even when tension vanishes on the time-slice. 

V. ORDER OUT OF DISORDER IN THE 
SYSTEM OF STRIPES. 

The main conclusion of the previous section is that a 
single quantum stripe, as defined through the sublattice 
parity order, is at best a very mildly fluctuating object. 
Considering a system of these Gaussian strings, at the 
moment one adds any interaction it has to be that long 
range order sets in. Algebraic order (of the single stripe) 
changes in true order in the presence of any perturbation, 
regardless its strength. Only recently exceptions have 
been identified (the quantum smectic, or gliding phase, 
Kivelson, Fradkin and Emery, 1998; Emery et at, 2000). 
However, these are only realized under specific circum- 
stances which are not found in the present context. 

Hence, any direct interaction between the stripes suf- 
fices to cause translational symmetry breaking in the sys- 



tem of stripes. There is no doubt that the stripes are 
interacting. They are charged and therefore they should 
exert Coulomb forces. In addition there are the Casimir- 
type forces in the spin system (Pryadko, Kivelson and 
Hone, 1998), as well as the elastic forces mediated by the 
lattice. 

Although mostly of academic interest, hard-core inter- 
actions (or non-intersection conditions) are special (Za- 
anen, 2000). These interactions are highly singular and 
a priori one cannot be sure that the hard-core interac- 
tion will play the same role as finite range interactions. 
Although there is no real good reason, it is appealing 
to assume that the stripe-stripe interaction contains a 
hard-core piece. One might want to be interested specif- 
ically in the question to what extent can stripes be a one 
dimensional sub-reality in two dimensional space. For 
this purpose alone one would like to keep stripes from 
intersecting. In addition, if one just wants to general- 
ize Ogata-Shiba to one higher dimension, one also better 
keep their hard cores attached to the charges. 

As discussed in section II, the charge sector of the 
Luttinger liquid of a strongly coupled Hubbard model 
is described by a hard-core bose gas. A most literal gen- 
eralization of this Luttinger liquid to 2+lD can be ob- 
tained by just a decompactification, in the same sense 
as used in fundamental string theory. In path-integral 
language, the hard-core bose gas corresponds with me- 
andering elastic worldlines, directed along the time di- 
mension, which cannot intersect (hard-core condition). 
At distances large compared to the lattice constant a 
one cannot see the difference between this system and 
a system of strings characterized by one more space di- 
mension which is curled up in a circle with compacti- 
fication radius i? ~ a, with the string wrapped around 
this extra dimension. Decompactification means that the 
compactification radius R ^ oo. What happens? The 
tiny string cylinders spread out in 2D worldsheets, cor- 
responding with elastic membranes, spanning the extra 
space dimension. The hard-core condition means that 
these worldsheets cannot intersect. This entity was called 
the directed string gas in 2-l-lD. The emphasis should be 
on directed because this decompactification construction 
gives rise to a constraint which is a-priori not completely 
general. On the time-slice the strings are directed along 
the extra dimension. A string starting at —oo in this 
direction always ends up at +oo in the same direction. 

It is a fundamental requirement of non-relativistic 
quantum- mechanics that worldlines/worldsheets are di- 
rected along the time direction. However, no general 
constraint of this kind acts in space directions, and in 
principle 'overhangs' or 'dislocations' (Fig. |lj), where a 
string for example starts out at —oo to end at — oo (for 
open boundaries), are in principle possible. 

A first objection could be that a single string can be 
subjected to directedness symmetry breaking, the sur- 
prise in the previous section. If the constituents are 
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already directed, the system will be definitely directed. 
However, although it is demonstrated that lattice strings 
of the previous section can acquire spontaneously a di- 
rection, there is no theorem available stating that lattice 
strings are always directed. Hence, one cannot claim that 
lattice string gasses are universally directed. 

However, there is an elegant argument available 
demonstrating that directedness is an unavoidable conse- 
quence of the dynamics in the system of hard-core elastic 
strings. This goes hand in hand with the demonstration 
that the directed string gas has to solidify (to break trans- 
lation symmetry) always. Exceptions are not possible. 
Hence, together with the physics discussed in the previ- 
ous section, the conclusion is that if the Ogata-Shiba ge- 
ometric squeeze prescription applies literally, long range 
order is unavoidable at zero temperature in 2+1 dimen- 
sions! 

Let us discuss the string-gas in more detail. The the- 
oretical problem is that due to the absence of a second 
quantization formalism the canonical methods of quan- 
tum mechanics are of no use for string problems. Hence, 
all what remains is the path integral formalism and in 
this formalism the string-gas problem corresponds with 
the statistical physics problem of elastic membranes em- 
bedded in 3D space subjected to a non-intersection con- 
dition, with the added constraint that the membranes are 
directed along one (imaginary time) direction. 

Let us step back, to reconsider the (seemingly) easier 
'compactified' version corresponding with directed, non- 
intersecting elastic lines in 2D. This is equivalent to the 
1-l-lD hard-core bose gas and it is well known that this is 
in turn equivalent to the problem of non-interacting spin- 
less fermions in 1+lD. This is of course a trivial problem 
and the freshman can calculate the density-density cor- 
relator of the fermion gas to find, 

nrnO =--— + — 13 

and the textbook will stress that these are the famous 
Friedel-oscillations, characteristic for any fermi-gas in 
any dimension. 

However, much later one learns that the spinless- 
fermion gas is just a Luttinger liquid characterized by 
a charge stiffness Kp = 2. In turn, since the observa- 
tions by Haldane (Haldane, 1981) and others it is clear 
that Eq. (|lj) has to do with algebraic long range or- 
der. Hence, the bosons order in a 1-l-lD crystal. This 
crystal is carrying phonons and the admixture of these 
phonons in the ground state change the true long range 
order in the algebraic order signaled by Eq. ([T^). This 
appears as a paradox: the Fermi-gas is mere kinetic en- 
ergy and how can this gas ever be a crystal? The reso- 
lution is that Fermi-statistics codes for a hard-core con- 
dition in the Bose language, and the hard-cores cause 
microscopic kinetic energy to become potential energy at 



large distances, driving the order. An interesting order- 
out-of-disorder mechanism is hidden behind the simple 
non-interacting fermions! 

This mechanism is well known in the statistical physics, 
addressing the problem of classical incommensurate flu- 
ids (domain wall fluids) in 2D (Pokrovsky and Talapov, 
1979; Coppersmith et ai, 1982). The argument goes back 
to work by Helfrich, 1978, actually on extrinsic curvature 
membranes in 2-l-lD, and was apparently reinvented in 
the community working on 2D incommensurate fluids. In 
the 1-t-lD context one can either use an intuitive argu- 
ment or a more rigorous self-consistent phonon method 
invented by Helfrich. In 1-l-lD one arrives at the same 
answer (at least qualitatively) but this is different for 
elastic strings in 2+lD, where the intuitive argument is 
flawed. Nevertheless, the intuitive argument is instruc- 
tive because it sheds light on the basic physics at work. 

This arguments is as follows for the 1+lD case. The 
hard-core bose gas at zero temperature corresponds with 
the statistical physics problem of a gas of non- intersecting 
elastic lines embedded in 2D space-time, which are di- 
rected along the time direction. The space-like displace- 
ment of the i-th worldline is parametrized in terms of a 
fleld 0i(T) (r is imaginary time) and the partition func- 
tion is {AI is the mass of the particle), 

drY^^idA^?. (14) 

supplemented by the avoidance condition, 

01 < (/)2 < ... < (l)N- (15) 

The hard-core condition Eq.(p^) renders this to be a 
highly non-trivial problem. 

At short distances the worldlines can freely meander. 
However, after some characteristic time-like distance, the 
worldlines will collide. In the statistical physics analogy, 
every collision costs an entropy ^ ks because the lines 
cannot intersect. Hence, these collisions raise the free en- 
ergy of the system and this characteristic free energy cost 
Ai^co/; ^ kBTucoU.- The density of collisions ricoii. is eas- 
ily calculated for the elastic worldlines. It follows from 
equipartitioning that the mean-square transversal dis- 
placement as function of (time-like) arclength increases 
like {[4>{t) — (/'(O)] ) = {h/M)T. The characteristic time 
Tc it takes for one collision to occur is obtained by im- 
posing that this quantity becomes of order <P where d is 
the average worldline separation, while the particle den- 
sity n ^ 1/d. A characteristic collision energy scale is 
obtained Ep ~ fi/T^ ~ {h'^/M)n'^. Ep is of course the 
Fermi-energy: it is the scale separating a regime where 
worldlines are effectively isolated [E > Ep, free parti- 
cles) from the one dominated by the collisions {E < Ep, 
Luttinger liquid). 
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At the same time, the entropy /kinetic energy cost gives 
rise to an effective repulsion between the world-lines, and 
this repulsion is in turn responsible for the ordering ten- 
dency. At large distances the precise origin of the repul- 
sion does not matter and one can simply assume that 
the entropic repulsion is like a harmonic spring and the 
spring constant can be estimated by taking the ratio of 
the characteristic energy (Ep) and the characteristic dis- 
tance d. In this way one finds a 'induced modulus' B 
associated with the compression of the hard-core 1+lD 
quantum gas, 



Bo - Ep/d 



Md^' 



(16) 



Asserting that at long wavelength the gas is described by 
the elasticity theory of a ID quantum crystal with spatial 
modulus Bq and mass density p ~ M /d, 

Seff^lJdTjdx [pidri^f + Boid^i^f] , (17) 

one recovers the spinless-fermion results, modulo prefac- 
tors of order unity. 

The more rigorous argument by Helfrich, 1978, starts 
out by assuming that the Bose-gas is described by the 
long wavelength action Eq. (p^). In the absence of the 
hard-core interaction Bq would be zero by definition and 
the free energy increases for a finite Bq because the fiuc- 
tuations are suppressed. Define a 'free-energy of mem- 
brane joining' as AF = F{Bo) - F{Bo = 0). At the 
same time, by general principle it has to be that the true 
modulus in the space direction B should satisfy {V is the 
volume). 



B^d 



,,d HAF{Bo)/V) 
dd^ 



(18) 



In case of the steric interactions, the only source of 
long wavelength rigidity is the fluctuation contribution 
AF. Therefore B = Bg and B can be determined self- 
consistently from the differential equation Eq. (|l8|). This 
method is not exact, because mode couplings are ne- 
glected. However, these mode couplings are important 
at short distances and they are therefore not expected to 
change the outcomes qualitatively. The ultraviolet only 
enters the answers through the short distance cut-off in 
the integrals, Xmm = V'i ^^'^ it appears that all the effects 
of these interactions can be absorbed in the fudge factor 
rj. Evaluating matters for the hard-core Bose gas, it turns 
out that it reproduces exactly the spinless fermion results 
ff 7? = V6 (Zaanen, 2000). 

The conclusion is that the algebraic translational or- 
der hidden in the hard-core Bose gas/spinless fermion 
problem in 1-l-lD can be understood as an order-out-of- 
disorder phenomenon in the equivalent statistical physics 



problem, which can be handled rather accurately, using a 
simple statistical physics method. The advantage is that 
the Helfrich method applies equally well to the string gas 
problem in 2-l-lD. In fact, it works even better! 

Let us first consider the directed string gas. The bare 
action of this string gas in Euclidean space-time describes 
a sequentially ordered stack of elastic membranes. Ori- 
enting the worldsheets in the y, r planes, the action be- 
comes in terms of the dispacement fields 4>i{y, r) describ- 
ing the motion of the strings in the x direction. 



Z^Ul.Uy^r I dMy,r)e 



S = J drdy ^ 



(19) 



again supplemented by the avoidance condition Eq. (|1E_ 
In Eq. (|l|), Pc is the mass d ensity and Sc the string 
tension, such that c = ^Sc/pc is the velocity. 

Let us now consider the intuitive collision-argument 
for this string gas. The mean-square transversal displace- 
ment now depends logarithmically on the worldsheet area 
A: ((A(/)(A))2) = h/{pc)\n{A). Demanding this to be 
equal to d^, the degeneracy scale follows immediately. 
The characteristic worldsheet area Ac for which on av- 
erage one collision occurs is given by h/{pc) ln{Ac) ~ d^ 
where Ac — c^T^ja^ in terms of the collision time Tc- It 
follows that Tc — {a/c)e^/'^^ and the 'Fermi energy' of the 
string gas is of order Ef^ = Ti/tc ~ (?ic/a) exp (— 1/2/i) 
where p is the coupling constant ('dimensionless K) of the 
string gas (Zaanen, Horbach and van Saarloos, 1996), 



pcd"^ 



(20) 



For a continuum description to make sense, p < 1 
and this suggests that the Fermi energy is exponentially 
small. However, it is finite and this is all what matters 
as we will see. 

One could be tempted to estimate the induced modulus 
by asserting B ^ Ep^ . However, contrary to the Bose- 
gas the above intuitive argument is qualitatively flawed 
and the reason is that the kinetic repulsions are no longer 
driven by the physics at the collision length scale. A sin- 
gle string is itself a 1-flD elastic entity, characterized 
by long wavelength fluctuations which are dangerous in 
the sense that these are responsible for changing true 
long range order in algebraic long range order. Helfrich's 
self-consistent phonon method shows that these long- 
wavelength single string fluctuations are also the ones re- 
sponsible for the induced modulus (Zaanen, 2000). Car- 
rying out the integrations one finds for the free-energy of 
membrane joining. 



AF 
~V~ 



irhc .B 5 
24^c d^ i 



In 



v'T.c d 
a2 B 



) + OiX'). (21) 
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expanding matters in the small parameter A = 
{'•/Ba)/{^\/Y.'qd). Since B is tending to zero, the loga- 
rithm is dominating and this term originates in the small 
momentum cut-off (long wavelength limit) in the integra- 
tion of the on-string fluctuations. 

The differential equation obtained by inserting Eq. 
( pl| ) in the self-consistency condition Eq. ( p^ can be 
solved and this yields, 



B = Ad^ 



(22) 



where A is an integration constant while fj, is the cou- 
pling constant defined in Eq. (po|). Hence, instead of the 
exponential of the 'naive' argument, a stretched expo- 
nential is found and this difference is entirely due to the 
logarithm in Eq. (^l|), finding its origin in the long wave- 
length on-string fluctuations. Hence, it is in this sense 
that the solidiflcation of the string gas is driven by the 
longest wavelength string fluctuations. 

Although the induced modulus is larger than naively 
expected, from a more practical viewpoint it is still quite 
small and it tends to be overwhelmed by the effects of 
finite range interactions. This reflects of course the fact 
that strings fluctuate much less than particles. However, 
we set out to demonstrate that long range order cannot 
be avoided in the string gas and for this purpose all what 
matters is that the modulus B is finite at zero temper- 
ature. This is a sufficient condition to exclude a zero- 
temperature proliferation of dislocations. In the absence 
of the dislocations (Fig. |lj) the string gas is sponta- 
neously directed and the directed gas solidifies always, 
as we showed in the previous paragraphs. 

The argument that dislocations cannot proliferate at 
zero temperature is quite nontrivial (Pokrovsky and Ta- 
lapov, 1979; Coppersmith et al. 1982). The string-gas 
theory Eq. (|l^) is generalized to finite temperature by 
compactifying the imaginary time axis with radius Rt = 
h/ksT. The non-proliferation theorem follows directly 
from the well-known result that a Kostcrlitz-Thouless 
transition (dislocation unbinding driven by thermal fluc- 
tuations) happens in this classical string gas at a finite 
temperature as long as the zero-temperature modulus is 
finite. Hence, dislocations are already bound at a finite 
temperature and they remain to be bound at zero tem- 
perature. 

A detailed analysis of the finite temperature case will 
be presented elsewhere (Mukhin, van Saarloos and Za- 
anen, 2001). The bottomlinc is that at finite tempera- 
tures one can simply use the high temperature limit of 
Eq. ( p^ (without the time direction), adding however 
the induced zero-temperature modulus ~ B{(j>i — </)j_(.i)^. 
This is nothing else than again the hard-core bose gas 
but now in its classical interpretation of thermally fluc- 
tuating elastic lines. The qualitative difference with the 
quantum case is that there is no longer a directedness 
constraint on the lines and in this classical gas disloca- 



tions can occur. If B = the remarkable result is that at 
any finite temperature dislocations are proliferated, while 
at the same time for any finite B the Kosterlitz- Thouless 
temperature occurs at a finite temperature Tkt ^ B. 

This has been discussed elsewhere at great length 
(Coppersmith et al., 1982) and let us just repeat the 
essence of the argument. The dislocations interact with 
long range, logarithmic forces which are set by the elastic 
moduli of the medium and therefore the energy associ- 
ated with free dislocations is logarithmic in the system 
size. At the same time, the entropy associated with free 
dislocations is also logarithmic and balancing these two 
yields the Kostcrlitz-Thouless criterion for the stability 
of the algebraic order. 



ady/Bx^c 
2^ 



> 1 



(23) 



Using the transfer-matrix the induced modulus Bt of 
the classical problem can be calculated exactly. Modulo 
prefactors this is Eq. ( pj| ) expressed in classical units. 
One finds \JBt^c ~ T which means that either the KT 
criterion is never satisfied (meaning that dislocations are 
always bound) or that the KT criterion is always satis- 
fied so that dislocations are proliferated at all tempera- 
tures. It turns out that the prefactors conspire in such 
a way that for two flavors of domains (our case) the sec- 
ond possibility is realized. This means that at any finite 
temperature dislocations always proliferate but they do 
so in the most marginal way. The entropic interactions 
driven by the finiteness of temperature are on the verge 
of beating the entropy of the dislocations but the former 
just loose. Any interaction other than this entropic in- 
teraction (including the quantum 'entropic' interaction) 
can tip the balance (Mukhin, van Saarloos and Zaanen, 
2001). Hence, adding a finite zero temperature B causes 
the Kostcrlitz-Thouless temperature to happen at a finite 
temperature. 

The conclusion is, remarkably, that nothing can keep 
the string gas away from solidifying at zero temperature 
(Zaanen, 2000). 



VI. THE QUANTUM MAGNETISM OF STATIC 
STRIPES. 

The magnetism of the stripe phase is relatively easy to 
study experimentally, and for this reason it is a relatively 
well developed subject. To put the remainder of this sec- 
tion in an appropriate perspective let us therefore start 
out with a sketch of the present empirical picture. 

It is a rather significant empirical fact that despite 
a high hole density well-developed antiferromagnetic or- 
der can be realized in doped cuprates (Tranquada et at, 
1995; Klauss et ai, 2000). It is well understood that 
a single hole is a strongly frustrating influence in the 
quantum-antiferromagnet. Since the spin system itself is 
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quite quantum-mechanical {S = 1/2) these frustrations 
give rise to the formation of a droplet of quantum spin liq- 
uid surrounding the hole (Dagotto, 1994). If these holes 
would stay independent, antiferromagnetic order would 
disappear at a very low doping. The very fact that Neel 
order has been demonstrated to persist in some systems 
to dopings as large as 20 % (Klauss et ai, 2000) should be 
taken as the leading evidence for the hypothesis of Sec- 
tion 11. Of course, it is also experimental fact that this 
big Neel order occurs when the charges organize in the 
stripes. However, in doing so the spin system becomes 
unfrustrated and this should be understood as the mani- 
festation of the Ogata-Shiba squeezing principle at work 
in 2-f ID. 

However, on closer inspection one finds that the stripe- 
antifcrromagnet is a more quantum-mechanical entity 
than the antiferromagnet of the half- filled insulator. 
Both NMR measurements (Hunt et ai, 1999; Curro et 
ai, 2000; Teitelbaum et ai, 2001) and neutron scattering 
(Tranquada, Ichikawa and Uchida, 1999) indicate that 
the spin-stiffness is smaller than the one at half-filling. 
It has been claimed that this should be due to a dilution 
effect: the exchange bonds connecting spins on opposite 
sides of the stripes (J') would be very small as compared 
to the exchange interactions inside the magnetic domains 
(J) which are in turn believed to be of the same magni- 
tude as the exchange interactions at half-filling. How- 
ever, for several reasons this cannot be quite the case. 
First, J' sets the scale for the overall incommensurate 
behavior and at energies larger than J' incommensurate 
spin fluctuations cannot exist (Zaanen and van Saarloos, 
1997). These fluctuations have been seen up to energies 
of - 40mey (Aepph et ai, 1997; Mook et ai, 2000) and 
this sets a lower bound to the value of J'. More directly, 
some inelastic neutron scattering data are available for 
the spin waves in a static stripe phase and these demon- 
strate that although the stiffness is strongly reduced the 
spin wave velocity stays large (Tranquada, Ichikawa and 
Uchida, 1999). This behavior is characteristic for the 
generic long wavelength physics of a Neel state which is 
on the verge of undergoing a quantum phase transition 
into a quantum-disordered state (Sachdev, 1999, 2000). 

The above observations are associated with the 214 
system. Recently, Mook et al, 2001, reported evidence for 
static stripes in the strongly underdoped 123 cuprates. 
However, they also claimed that although charge order 
is established, the spin system is apparently quantum 
disordered. The incommensurate spin fluctuations are 
seen only above a small but finite (— 3 meV) energy. 
This is not surprising. It is well understood that the bi- 
layer couplings as they occur in 123 are a factor promot- 
ing quantum spin fluctuations (Millis and Monien, 1993; 
van Duin and Zaanen, 1997). Since the spin system in 
the single layer 214 cuprates is already on the verge of 
quantum-melting, these bilayer couplings could easily tip 
the balance. 



NMR measurements have shown that the actual 
asymptotic spin-ordering process is highly anomalous 
(Hunt et ai, 1999; Curro et al, 2000; Teitelbaum et al, 
2001). It appears that slow spin fluctuations (MHz scale) 
show up at the temperature where the scattering exper- 
iments indicate a freezing behavior (— 70 K), to con- 
tinue down to the lowest measured temperatures (400 
mK). These fluctuations are at present not at all un- 
derstood. However, although the case is definitely not 
closed, it appears that the spin dynamics on a larger en- 
ergy scale fits quite well the expectations of the generic 
field theory describing the long wavelength dynamics of a 
coUinear quantum-antiferromagnet close to its quantum 
phase transition. All what matters is the symmetry of the 
order parameter (0(3)) and the dimensionality of space- 
time: this generic theory is the 0(3) quantum non-linear 
Sigma model in 2-hl D (QNLS). 

Several excellent treatises are available, both on the 
introductory (Sachdev, 2000) and the advanced level 
(Chakravarty, Nelson and Halperin, 1989; Sachdev, 
1999), on the physics near quantum phase transitions. 
Let us therefore limit ourselves to the bare essence. 
It is well understood that the non-frustrated Heisen- 
berg quantum-antiferromagnet defined on a bipartite 
lattice does not suffer from Marshall sign problems. 
Stronger, the long wavelength dynamics in the semi- 
classical regime is free of Berry-phases and it can there- 
fore be described with the simple QNLS (Fradkin, 1991), 

Z = J Vn6{\n\ - 1)6"^ 

S= — [ d'x f dT{{drnf + (Vn)2) (24) 
9o J Jo 

in scaled variables, such that the spin wave velocity is 
one. n is a three component vector of fixed length and 
Eq. (|2j) is nothing else than the theory of a classical 
Heisenberg spin system embedded in 2-t-l dimensional 
Euclidean space time (Chakravarty, Nelson and Halperin, 
1989). At zero-temperature (/? oo) this becomes pre- 
cisely equivalent to the classical Heisenberg problem in 
3D. Hence, for small bare coupling go (low temperature 
in the classical problem) Neel order is established. At 
a critical value gj? a second order phase transition oc- 
curs to a strong coupling, quantum disordered state. In 
the classical problem the spin correlators decay exponen- 
tially in the disordered state and this means that the real 
time dynamics of the quantum problem is characterized 
by a dynamical mass gap in the mode spectrum. Right 
at the critical point the dynamics is scale invariant while 
in the proximity of this point the same is true up to a 
length/time scale (Josephson correlation length) where 
the system finds out that it either gets attracted to the 
Neel fixed point or to the disordered state. At higher 
energies and larger momenta the system still behaves as 
if it is at its critical point, until a crossover is reached 
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below which one sees the dynamics associated with the 
Neel state (zero-modes) or the disordered state (massive 
triplet excitations). The behavior at finite temperature 
is especially interesting. Finite temperature means com- 
pactification of the imaginary time axis with compactifi- 
cation radius Ti/ksT. In the critical regime this breaks 
scale invariance meaning that the dynamics is character- 
ized by a characteristic time and length Tc ~ 2c ~ h/ksT. 

The question arises what this means for the real time 
dynamics - the Wick rotation is a remarkably counter- 
intuitive affair. The answer is Sachdev's achievement 
(Sachdev, 1999). At times r < Tc the quantum field 
theory itself generates through the analytic continua- 
tion a classical relaxational dynamics characterized by 
the 'quantum limit of dissipation': this relaxation time 
is as short as it can be, namely h/ksT. At shorter times 
the zero-temperature critical dynamics is recovered, char- 
acterized by a cusp-like dynamical susceptibility reflect- 
ing the criticality in space-time. Away from the critical 
point, it is the same at high temperature, high frequency 
and large momenta, but at temperatures low compared 
to the Josephson scale the mode-excitations associated 
with the stable fixed points dominate the low-frequency 
end of the spectra. 

This theoretical picture seems to explain the most 
salient features of the spin dynamics in the cuprates, 
at least if one identifies this spin dynamics with the 
incommensurate spin phenomena which are most nat- 
urally interpreted as being related to the stripe Neel 
state. The spin dynamics of the static stripes would then 
be interpreted as reflecting the classical sector, that of 
the fully developed superconductors at low temperatures 
with that of the disordered massive state (characterized 
by a spin gap in the incommensurate spectrum), while 
the normal state well above the ordering temperatures 
would be related to the critical regime. A crucial assump- 
tion is of course that some form of spin-charge separation 
takes place and this is in fact clearly excluded by the ex- 
periments: the spin gap appears at the superconducting 
transition. Hence, this interpretation is at best only part 
of the explanation. 

However, in the context of static stripes the situation 
is more clear. The interpretation in terms of the stripe- 
antiferromagnet being on the verge of undergoing the 
quantum-phase transition rests on data obtained below 
the charge ordering temperature, and in this context a 
spin-only interpretation is more reasonable, while it adds 
credibility to the notion that QNLS has something to do 
with high Tc. 

The next question is, what is the source of the strong 
quantum fluctuations in the stripe antiferromagnet? This 
can have several reasons. A first obvious possibility is, 
in the language of this paper, 'local violations of the 
Ogata-Shiba squeeze' - by longer excursions away from 
the stripes holes can cause local violations of the bipar- 
titeness of the squeezed lattice. These correspond with 



local spin-frustrations and they should therefore enhance 
the collective fluctuations seen at long wavelength. How- 
ever, even in the case that the squeezed lattice is per- 
fectly bipartite it appears still possible to end up with a 
quantum disordered stripe anti-ferromagnet. After this 
detour we are back at the 'mean-fleld theory of sublattice 
bipartitencss'. 

There is no reason to assume that the strength of the 
exchange bonds in squeezed space is the same every- 
where. The exchange interaction (J') between the spins 
on opposite sides of the stripe is caused by a microscopic 
dynamics (hole motions) of an entirely different kind than 
the superexchange which is responsible for the spin-spin 
interactions inside the magnetic domains (J). Although 
nothing is known for certain, it is generally expected that 
J' < J. If J' would vanish, the 2+ ID spin system associ- 
ated with the ordered stripes would be cut into indepen- 
dent 1-l-lD spin ladders with an effective width set by the 
stripe separation and the details of the stripe ordering. 
If these ladders have an even width, it is well understood 
(Dagotto and Rice, 1996) that for J' is zero the system 
would be characterized by a spin gap. This spin gap of- 
fers protection for the quantum disordered ladder state 
at flnite J': a critical value of has to be exceeded 
before classical Neel order can emerge. For domains of 
uneven width, and the case of uneven spin, the spin sys- 
tem on every ladder is a Luttinger liquid (Dagotto and 
Rice, 1996) and any J' will suffice to cause long range 
order. 

This ladder notion acquires an additional significance 
in the light of recent theoretical works addressing the 
microscopic mechanism of stripe formation. Much of 
the earlier work was based, implicitly or explicitly, on 
the large S limit. Here the quantum spin fluctuations 
are neglected completely and one finds the site ordered 
stripes as they first appeared in the Hartree-Fock calcula- 
tions by Zaanen and Gunnarsson (Zaanen and Gimnars- 
son, 1989), and variations thereof. Recently, Voita et al. 
(Voita and Sachdev, 1999; Voita, Zhang and Sachdev, 
2000) considered a limit which is in a sense opposite 
to large S: t — J type models are characterized by a 
global SU{2) symmetry and these can be generalized to 
a Spl{2N) symmetry. By sending N — > oo, keeping 5* fi- 
nite, saddle points can be identified characterized by ex- 
ceedingly strong quantum-spin fluctuations when viewed 
from the large S side. In this limit, the 'spin' system 
is generically unstable towards the formation of spin- 
Peierls (or 'valence bond') phases (Read and Sachdev, 
1989; Sachdev, 1999, 2000). Nearest-neighbor spins form 
pair wise singlets and these singlets are stacked in ladder- 
like patterns on the 2D planes. Voita et al. showed that 
the large A'^ saddle-points also correspond with stripe 
phases as long as the hole density is not too large. The 
magnetic domains appear as Peierls-ordered even-leg lad- 
ders while the stripes are like highly doped two-leg lad- 
ders ('bond ordered stripes'). The additional benefit is 
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that at large A'', these charge stripes are generically su- 
perconducting while uniform d-wave superconductivity 
takes over at large dopings. 

The most trustworthy microscopic calculations avail- 
able at present are the numerical DMRG studies by 
White and Scalapino (White and Scalapino, 1998). 
These calculations indicate that the stripes of the t — J 
model are somewhere in the middle of large S and large 
N . On the one hand, the stripes are bond centered 
and a case is made on basis of the numerics that these 
stripes have a tendency to become superconducting. On 
the other hand, diagonal site centered filled stripes are 
nearby in energy and these are quite like the Hartree- 
Fock stripes. In addition, it appears to be easy to stabi- 
lize Neel order and this order is characterized by a strong 
anti-phase boundarieness. 

To obtain a better understanding of the long wave- 
length quantum magnetism of the stripes we studied our- 
selves in great detail ladder-like spin models numerically 
using quantum Monte-Carlo (Tworzydlo et al., 1999). 
These can also be considered as being representative for 
the generic spin-dynamics of the two dimensional system 
where the Ogata-Shiba squeezing applies literally while 
the charge stripes are static. Define a S ~ 1/2 nearest- 
neighbor Heisenberg Hamiltonian on the 'squeezed' bi- 
partite lattice. On this lattice the stripes corresponds 
with a regular array in, say, the x direction of lines ex- 
tending along the y-direction which are centered on the 
links of the lattice. The exchange interaction is J every- 
where, except for the links which are 'cut' by the stripes 
where the exchange interaction is J'. Besides temper- 
ature, the free parameters are (a) a = J' /J, the ratio 
of the stripe-mediated exchange interaction and the su- 
perexchange, and (b) the number of sites Niegs separating 
the J' links in the y direction: the stripe separation. 

This is a simple bond-dilution Heisenberg model which 
can be studied to any desired accuracy using the novel 
cluster-loop algorithm quantum Monte-Carlo method. 
Obviously, for a > 1 Neel order cannot be avoided at 
zero-temperature and the interest is in what happens for 
small a. Let us therefore first consider a ^ 0. The spin 
system is qualitatively different in this limit for even and 
uneven Niegs ■ For uneven Niegs it corresponds with a dis- 
connected system of spin-ladders with an uneven number 
of legs and it is well known that these ladders renormalize 
in Luttinger liquids. The ladder-to-ladder coupling J' is 
in this situation always relevant and for any finite a the 
ground state will exhibit long range order. The way this 
ground state is approached as function of temperature 
is remarkably simple and is illustrated in Fig. (^) for 
the one-leg and three-leg cases which appear to behave 
in a near-identical way. At finite temperature all correla- 
tions arc short ranged but several correlation lengths and 
associated characteristic temperatures can be identified. 
First, as function of decreasing temperature the correla- 
tion length ^1 in the ladder direction will rapidly increase 
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FIG. 6. Crossover temperatures in units of J (see text) as 
a function of the anisotropy a for the coupled one- and three 
leg (a), as well as the two leg (b) spin- ladder models. The 
lines and points refer to the analytical- and numerical results, 
respectively, for the various scales. 

and it is for small a by far the largest length. At distances 
less than the spin-correlators in this direction will ex- 
hibit the algebraic correlations of the zero temperature 
case. For small but finite a there will be a temperature 
Tq where the ID correlations become so strong that even 
a small a will sufRce to cause the spins to correlate in 
a 2D fashion. This dimensional crossover temperature is 
defined by the temperature where the correlation-length 
in the x direction becomes of order of twice the width of 
the ladder. We found that Tq ^ a for small a which ap- 
pears to be consistent with the scaling theory of Affleck 
and Halperin, 1996. 

At lower temperatures the spin-system is 2-1-1 dimen- 
sional and its long wavelength dynamics should be char- 
acterized by the universal behaviors which follow from 
the 0(3) QNLS model. At a given temperature T < To 
there are three possibilities: 

(a) the spin correlation length tends to saturate to a fi- 
nite value at zero-temperature. This signals that the bare 
coupling constant go > Qq-, and the system flows to the 
quantum disordered zero temperature state. 

(b) The correlation length behaves in the renormalized 
classical fashion, meaning that the spin system is un- 
dergoing only thermal fluctuations although the spin- 
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stiffness is smaller than expected because of the influence 
of the quantum fluctuations at shorter scales. In this 
case one expects for the temperature dependence of the 
correlation length ^(T) = exp(T*/r)/(2T* + T) where 
T* = 2'!rps in terms of the renormalized stiffness ps ■ 
(c) The temperature is larger than either the zero tem- 
perature spin-gap associated with case (a) or the spin- 
stifFness of case (b) so that the 2-l-lD spin system still 
resides in the quantum critical regime. In this case 

^ 1/T, revealing that temperature sets the scale. 

As it turns out, in the uneven- leg cases the system 
jumps directly into the renormalized classical regime at 
the moment it finds out that it it becomes two dimen- 
sional. Hence, at higher temperatures the physics is that 
of decoupled 1-l-lD subsystems and these discover that 
they are on the way to a classical 2D Neel state at the mo- 
ment that the temperature is low enough such that these 
1+lD subsystems start to correlate in a 2 -I- ID fashion. 
It is even so that the renormalized spin stiffness T* ~ T*^, 
the dimensional crossover temperature. 

The even-leg case is far more interesting and this is 
illustrated in Fig.^ for the two-leg ladders. Even leg 
ladders in isolation (a = 0) are well known to exhibit 
a spin gap. This phenomenon is probably best under- 
stood as a consequence of the fact that the Peierls states 
are the natural competitors of Neel order for lattice- 
quantum antiferromagnets (Read and Sachdev, 1989; 
Sachdev, 1999, 2000). Define a Peierls order parameter 
field which amounts in the two leg ladder case to stronger 
exchange bonds along the rungs of the ladder and weaker 
bonds along the legs. Obviously, if this difference is large 
enough two spin singlets are formed along the rungs, and 
the overall state of the system can be viewed as a sim- 
ple row of these rung-singlets and the system has a spin 
gap. Upon reducing this difference to zero, the spin gap 
stays finite and therefore the state of the spin ladder with 
uniform exchange couplings is adiabatically connected to 
the Peierl's state where translational symmetry is explic- 
itly broken. For four leg-, six leg ladder, etcetera, cases 
the same argument can be used, except that one is now 
dealing with 2 and 3 parallel rows of singlets, respec- 
tively. The spin gap decreases rapidly with the number 
of legs in the absence of frustrations. However, adding 
frustrations it appears that a stable phase exists where 
this gap stays open even if the number of legs approaches 
infinity (Read and Sachdev, 1989; du Croo-de Jongh, van 
Leeuwen and van Saarloos, 2000). It is noted that this 
logic does not apply to uneven leg ladders because of the 
relevancy of the topological phases associated with the 
uneveness of the total spin on every rung. 

Let us now consider the cross-over diagram as func- 
tion of a and temperature for the coupled two-leg ladder 
problem. For a of order 1 the spin system is just like 
the uniform Heisenberg problem on a square lattice. At 
a temperature T/ J ~ 1 the correlation length is of order 
the lattice constant and from Fig. Hb one infers the well 



known result that this Heisenberg problem (of relevance 
to the half-filled case) is so far away from the critical 
point that before the crossover to the quantum critical 
point is reached the correlation length has already hit 
the lattice constant: at all temperatures this system is in 
the renormalized classical regime. This changes when a 
is reduced. As in the uneven leg cases, the dimensional 
cross-over temperature Tq decreases well. However, the 
renormalized stiffness of the 2D renormalized classical 
regime decreases more rapidly and a window opens up of 
genuine 2-1-1 dimensional quantum criticality. This is of 
course due to the quantum phase transition to the Peierl's 
state occurring at ac ~ 0.3. It is noticed that for the 
two leg ladder magnetic domains this transition occurs 
while J' is still quite substantial. In the proximity of the 
quantum-critical point one is dealing with the competi- 
tion between the Neel state and a truely 2-1-1 dimensional 
spin-Peierls instability which is just helped by a partial 
explicit breaking of translational symmetry. It is only 
for a < 0.2 that the one dimensional on-ladder spin dy- 
namics gets protected: the on-ladder correlation length 
never grows large enough to fulfill the conditions needed 
to cause a 2-f ID spin dynamics, and in this regime one 
is dealing with decoupled ladders. It is noticed that all 
these behaviors have been reproduced in quite some de- 
tail using the scaling theory for the spatially anisotropic 
0(3) quantum non-linear sigma model by van Duin and 
Zaanen, 1998, extending earlier work by Castro-Neto and 
Hone, 1996. 

VII. UNIFYING SPINS AND STRIPES: A 
SIMULATION. 

Starting out with the Ogata-Shiba principle of section 
II, we deduced a general problem and subsequently we 
studied several aspects of this problem separately. Of 
course, the whole can be more than the sum of its parts. 
What happens when the lattice strings of section IV form 
the interacting string system of section V, which is in 
turn communicating with the quantum-antiferromagnet 
of section VI? The answer is the phase diagram shown 
in Fig. [t] (Osman, 2000). The control parameters are 
the a of the previous section parametrizing the relative 
strength of the stripe-mediated spin-spin coupling and 
the hopping parameter t of section IV, controlling the 
strength of the stripe quantum fluctuations. Finally, the 
stripe density matters and the phase diagram is the one 
for a density such that the domains form two-leg ladders 
when the stripes are static. We also looked at differ- 
ent densities. The topology of the phase diagram is the 
same for all commensurate stripe densities such that the 
magnetic domains form even-leg ladders. For increasing 
ladder width the parameter regime where phase I is sta- 
ble shrinks rapidly, while it disappears completely for all 
other densities. In all other regards the resulting phase 
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FIG. 7. Phase diagram of the model describing quantum 
stripes fluctuating through a quantum anti-ferromagnet, for 
the particular case that stripe are separated at average by 
two sites, a is the spin-anisotropy parameters and t is the 
stripe kinetic energy, both measured in units of J. Phase I 
corresponds with a charge-ordered and spin quantum disor- 
dered stripes, while regions II and III correspond with full 
charge and antiferromagnetic order. However, regions II and 
III are separated by a sharp crossover associated with the 
single stripe unbinding transition. 



diagram is generic. 

There are actually only two phases: the charge-ordered 
incompressible stripe antiferromagnet found for small a 
and t (phase I) and a fully charge and spin ordered phase 
(II and III). However, although 'phases' II and III are in- 
distinguishable at long wavelength their short distance 
nature is very different. The line separating these two 
phases is no more than a cross-over line but it is a very 
sharp cross-over which was actually detected in the nu- 
merics using the Binder criterion (Binder, 1997) which is 
devised for phase transitions. This cross-over line corre- 
sponds with the single string unbinding transition of sec- 
tion IV. Hence, in phase II a single stripe breaks already 
translational symmetry and the 2D order in the local- 
ized stripe system is driven by stripe-stripe interactions 
mediated by the spin system. In phase HI every indi- 
vidual stripe is delocalized but the order-out-of-disorder 
mechanism as discussed in section V takes over, causing 
again translational symmetry breaking. However, this is 
a very weak order. Although the numerical calculations 
confirmed this mechanism it is in many cases very hard to 
detect. The novelty is that the interplay of the quantum 
spin dynamics of section VI and the stripe fluctuations 
gives rise to a next surprise. Sufficiently strong stripe- 
quantum fluctuations restore the Neel order even if the 
spin system of the static stripes would he quantum disor- 
dered. Although this Neel order is weak when measured 
by the standard spin-correlator, the hidden spin-order as 
measured by the topological correlator Eq. (H) as dis- 
cussed in section II approaches the magnitude of order 



as established in the Heisenberg system of the half-filled 
insulator, but only so when t is large. The conclusion is 
that by just imposing the Ogata-Shiba squeezing condi- 
tion, quantum fluctuations turn into agents, mediating 
every order which can be realized in this system. 

Let us discuss in more detail the model and the numer- 
ical simulations behind this phase diagram. The model 
is most easily constructed in squeezed space. Hence, we 
consider a two dimensional square lattice with a Heisen- 
berg 5 = 1/2 spin on every site. These spins interact via 
nearest- neighbor interactions J = 1 except along 1-1- ID 
connected trajectories defined on the links between the 
sites which are characterized by a weaker exchange in- 
teraction J' = a J, see Fig. (||). Different from the 
model in Section V these trajectories can now have arbi- 
trary shapes (see Fig. and by unsqueezing the lattice 
one recovers the lattice strings. Notice that the connect- 
edness constraint that the holes have to be nearest- or 
next-nearest-neighbors translates on the squeezed lattice 
in the constraint that the J' links have to be nearest- 
or next-nearest-neighbors on the link lattice. Notice also 
that for diagonal stripe configurations, nearest-neighbor 
spin bonds emerge on the squeezed lattice, correspond- 
ing with next-nearest-neighbor spins on the unsqueezed 
lattice (crosses in Fig. ||) and these particular squeezed 
lattice links carry therefore no spin-spin interaction. As 
we will discuss later in more detail, these missing ex- 
change interactions are equivalent to the parameter K in 
the string model of Section IV, expressing the energy dif- 
ference between nearest- and next-nearest-neighbor hole 
bonds. All other curvature parameters (the L's) are set 
to zero. 

The J' strings are quantized in the same way as the 
strings of section IV. The J' links can hop to nearest- 
neighbor positions on the link lattice as long as these 
moves do not violate the connectedness constraint. We 
learned in section HI that this combination always leads 
to the directedness symmetry breaking. Hence, undi- 
rected string configurations do not have to be taken into 
account, which simplifies the calculations considerably. 
At the same time this means that dislocations cannot 
occur which represents a serious limitation. However, 
this limitation follows from the 'mean-field' requirement 
that the squeezed lattice has to be partite everywhere. 

Finally, we impose a hard-condition by requiring that 
J' bonds cannot meet on the same link of the same lat- 
tice. As can be easily checked, this means in the un- 
squeezed space that two stripes should be separated by 
at least one spin site. The hard-core is therefore spread 
out and this does not matter, although one should be 
aware of it when the stripe density is high. 

This defines the model and let us for completeness for- 
mulate the Hamiltonian explicitly. Define coordinates 
(x,y) for the sites on the squeezed lattice. Consider the 
stripes to be directed along the y direction. Define hard- 
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FIG. 8. The model describing fluctuating stripes in the quantum antiferromagnet. Starting with various connected stripe 
configurations, one arrives at a spin-only model in squeezed space, characterized by fluctuating lines of modified exchange 
bonds. Notice that in squeezed space exchange bonds disappear (crosses) when the original stripe configuration contains 
diagonal pieces. 



core particles a^^^^^ = which live on 

the links of the squeezed lattice where (x, y) labels the 
link connecting the site {x,y) with the site (a; + 1 , y) . The 
Hamiltonian is, 



+J(1 - (1 - a)n(^x y))S(x,y) ■ S^x+i^y) 
+,/(! - {nt 

X ,y)'^{x — l,y-\-l) 
~^'^{x-l,y)'^{x,y+l))S(x,y) ' S(^x,y-\-l)) 



(25) 



The first term corresponds with the stripe kinetic en- 
ergy where V is the projection operator of Eq. (||), but 
now defined on the links, ensuring that hole hoppings 
do not break up the strings. The second term changes 
the exchange J into J' when the particle is present on 
the bond, and the last term takes care of the missing ex- 
change interactions associated with the diagonal stripes 
(the crosses in Fig. ^). 

We studied this model numerically using a 'hybrid' 
Monte-Carlo algorithm for the spins and the strings. The 
updating can be done quite independently in these sub- 
sectors because strings and spins live on different (site 
vs. link) lattices. The spin system was simulated us- 
ing the loop-cluster algorithm (Evertz, Lana and Marcu, 
1993) and improved estimators (Wiese and Ying, 1994) 
which were also used in the simulations discussed in sec- 
tion VI. The string system was simulated using the same 
world-line algorithm as employed for the single string 
problem, which is described at length in Eskes et ai, 
1998. Although this latter algorithm is a conventional 
Monte-Carlo algorithm, and therefore not as efficient as 
the loop-cluster algorithm used for the spin system, the 
dynamics is relatively simple and we could simulate sys- 
tems of up to 20 strings on a squeezed lattice of 60 x 60 
sites. To give some idea, in Fig. ^ a snap-shot is shown 
of a time-slice for a system of 12 strings which are at av- 
erage separated by four-leg ladder domains, represented 
in unsqueezed space. 

To determine the phase boundaries and cross-over lines 
in the phase diagram Fig. ^ more accurately we used 



the Binder parameter technique (Binder, 1997). This 
parameter is the reduced fourth order cummulant defined 



B = 1 



{O* 



3(02 



(26) 



where O is the order parameter of the ordered phase. 
This parameter behaves differently in the ordered and 
disordered state. In the ordered state (O^) — (0^)2 g^j^^j 
therefore B = 2/3 while in the disordered state char- 
acterized by Gaussian fluctuations (O*) — 3(0^)2 g^^^ 
therefore B — Q. This technique also turned out to be 
useful for finding the cross-over line between regimes II 
and III, indicating that the order parameter is suddenly 
strongly reduced at this cross-over. 

Let us now turn to the phase-diagram. The case t = 
just corresponds with the spin-only physics discussed in 
the previous section and the issue is what happens at 
finite t. A particular simple case is where a = 1. The 
spins and strings move independently except for the last 
term in Eq. (p5|). We already indicated that this term 
is equivalent to a finite /C parameter in the string-only 
problem. A positive and sufficiently large /C localizes a 
single string along the 'vertical direction' and only when 
t exceeds a critical value the string unbinds from the lat- 
tice. This single string unbinding transition is responsible 
for the sharp cross-over between the II and III regimes 
in the diagram Fig. |^. The critical t can be easily esti- 
mated. A kink in the string will break a spin-bond along 
the y-axis and this will cost an energy equal to the en- 
ergy per bond in the pure Heisenberg spin system. The 
energy per site has been calculated by many groups to be 
—0.6692 J. Given that for every site there are two bonds, 
and given that the cost of a kink K. is equal to that of a 
spin bond, we find K. = 0.34J. Eskes et al. found that 
the transition from the flat- to Gaussian string phases 
occurs when JC = O.Ttc and we therefore estimate the lo- 
cus of the cross-over to be at tc = 0.34 J/0. 7 = 0.49 J, in 
striking agreement with the simulations. 

For t < tc there is definitely a very strong stripe or- 
dering tendency. Every stripe is like a rigid rod, and any 
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FIG. 9. A snapshot of a timeslice of the quantum Monte-Carlo simulation for a system of 12 stripes which are at average 
four sites apart. The black lines are the stripes and the light- and middle gray shades (yellow and red in the color version) 
indicate opposite orientations of the hidden spin-order parameter (staggered spin on the squeezed lattice). 



stripe-stripe interaction will lead to translational symme- 
try breaking towards a periodic state. We will in a mo- 
ment discuss the nature of the stripe-stripe interactions 
mediated by the spin system. This is different in the 
regime t > tc- The stripes by themselves delocalize and 
the ordering tendency is now driven by the order-out-of- 
disorder effects discussed in section V. This is a very weak 
order and we actually did not manage to detect it in the 
simulations for densities where the stripes are separated 
by two-leg ladders at average. This in fact illustrates 
vividly the order-out-of-disorder mechanism. Because of 
the 'smeared out' hard-core of the stripes, the system is 
very dense at this stripe separation and the fluctuations 
are therefore strongly reduced. Intuitively one would ex- 
pect that there would be even less tendency towards order 
if the stripes are placed further apart, but the opposite 
happens. We calculated the charge-structure factor on 
the unsqueezed lattice for a system of 12 stripes which 



would have a periodicity in the ordered state of 5 lattice 
constants (domains are 4 legs wide, Fig. ||) at a large 
t = 8J. As can be seen from Fig. |l^, this structure 
factor is characterized by a quite sizable '2e' charge or- 
der peak, relative to the lattice Bragg peak located at 
the origin. Although this represents a striking qualita- 
tive confirmation of the string-gas order-out-of-disorder 
mechanism (Zaanen, 2000) we did not attempt to further 
quantify these matters because of the rather serious size 
limitations of our simulations. 

Let us now turn to the interesting aspects of the in- 
terplay between the spin- and the string system which 
emerges when a becomes small. A first issue is that the 
spin system mediates interactions between the stripes, 
even when the stripes are static (small t) . This has been 
studied in quite some detail by Pryadko, Kivelson and 
Hone, 1998. Assuming small a, the stripes correspond 
with reflecting boundaries for the spin waves and as a 
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other sources, like electron-phonon coupling and direct 
Coulomb interactions. 




FIG. 10. Plot of the charge-charge structural factor for a 
lattice of 48 x 48 at the point a — 0.05 and t — 8.0 inside 
phase III. The average distance between the strings is 4. 



result attractive Casimir forces arise which would render 
the stripe system unstable towards phase separation in 
the absence of compensating repulsive long range forces 
(like, e.g., direct Coulomb interactions). This Casimir 
potential falls off like V{d) ^ 1/d where d is the stripe 
separation. However, this long-wavelength analysis is not 
quite applicable to the most relevant cases where stripes 
are only a few lattice constants apart. In this situation, 
the ladder effects as discussed in the previous section are 
expected to become also quite important with regard to 
the stripe-stripe interactions. We studied this problem 
by inserting two static horizontal stripes separated by d 
spin sites on a large Heisenberg lattice, comparing its en- 
ergy with that of a single stripe, and the energy of the 
pure spin system. The stripe-stripe interaction energy 
per unit length of stripe is by definition. 



V^nt{d) = {E^id) + Eo{d) ~ 2E,id))/Ni 



(27) 



where £"„ is the total energy of the system with n stripes 
while Ni is the stripe length. We have calculated Vint 
for both a = 0.2 and a — 0.08 and the results are shown 
in Fig. 11. The major surprise is that even for stripe 
separations as small as two lattice constants these inter- 
actions are very weak, ~ 0.02J, to become even more 
minute at larger distances. The two leg ladder case ap- 
pears to be exceptionally stable, while the calculations 
suggest that even separations are always more stable 
than the uneven cases, as expected from the presence 
of a spin gap in the former. For instance, if all other 
forces could be neglected, these spin-mediated interac- 
tions would render a stripe system with d = 3 to be- 
come unstable to a stripe 'density wave' characterized by 

d = 2,4: 3-3-3-3 > 2-4-2-4 . 

However, it appears that in reality these quite feeble 
forces would be easily overwhelmed by interactions from 




FIG. 11. Induced interaction between two static flat 
stripes, d is the distance between the two stripes and e is 
the energy per unit length of stripe. 



We discussed in the previous section that static stripes 
separated by even leg ladders give rise to a quantum dis- 
ordered phase for small a. An interesting result is that 
strong stripe fluctuations restore the Neel state, and the 
gap of the incompressible spin phase (phase I) vanishes 
even for the two leg-ladder case before the single stripes 
unbind from the lattice (see Fig. In hindsight this 
is not so surprising. The 'hole' motions will cause kinks 
in the stripes and this means that the two leg-ladder 
is locally destroyed: over some length the spin domain 
becomes a 1 leg or 3 leg ladder (Fig. |l^). In this re- 
gard, it is interesting that the single stripe unbinding 
crossover increases substantially at small a's, more or 
less tracking the magnitude of the spin-gap in the in- 
compressible regime. Apparently this crossover scale is 
no longer completely due to the simple missing exchange 
bond mechanism we discussed earlier. We also learned 
that the energy associated with changing the stripe sep- 
aration locally is quite small (previous paragraph) and 
actually too small to explain the upturn of the cross-over 
line. Therefore, this upturn is caused by some non-trivial 
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FIG. 12. (a) Bound kink-antikink pairs in a coupled two-leg 
ladders stripe system, (b) When the kinks unbind strings of 
odd-leg ladders develop corresponding with a confining po- 
tential. It is believed that this mechanism is responsible for 
the upturn of the single string unbinding cross-over at small 
a. 

mechanism. We suspect that this is of the kind as illus- 
trated in Fig. [T^. In order for a single stripe to unbind, 
kinks should deconfine: isolated side steps of the strings 
should proliferate freely in the vacuum. From Fig. ^ 
one infers that for this too happen, strings of one-leg 
and three-leg ladders have to be created and these exert 
clearly a confining force on the kinks and such a confin- 
ing force is not considered in the simple string model of 
Eskes et al, 1996, 1998 . 

With regard to the recurrence of Neel order, this de- 
struction of two-leg ladderness is not all. In close analogy 
with the 1+lD case, at large t the holes are fluctuating 
rather freely and these hole motions induce a direct spin- 
spin interaction: let one hole hop back and forth and 
one directly infers that this causes a strong preference 
for the spins on both sides of the stripes to be antipar- 
allel. Hence, for t large as compared to J these hole 
motions average away the difference between J and J' 
and at large distances the exchange interactions can be 
taken to be uniform. Hence, A J = J — J' is in this sense 
an irrelevant operator and one recovers a notion of spin 
charge separation which is quite like the one encountered 
in the 1-l-lD Luttinger liquid. The difference is of course 
that in 2-1- ID the spin system orders. 

In squeezed space this charge-fluctuation induced Neel 
order is very easy to observe. The topological spin corre- 
lator, Eq. which measures this order in unsqueezed 
space (by 'dividing out' the anti-phase boundarieness at- 
tached to the charge) is equivalent to the direct spin cor- 
relator in squeezed space. At the same time, starting in 
squeezed space the direct spin correlator as measured in 
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FIG. 13. The behavior of the topological spin-spin corre- 
lator (pluses) and the direct staggered spin-spin correlator 
(crosses) as function of distance on a 24 x 24 squeezed lattice 
at a two- leg ladder stripe density {a = 0.05, t = 2.0). 



unsqueezed space is easily computed by reinserting the 
antiphase boundaries using a simple algorithm. In Fig. 

we show typical results for the direct and topologi- 
cal spin correlators in unsqueezed space in region HI of 
the phase diagram Fig. 0, calculated for a density cor- 
responding with two-leg spin ladders. It is seen that the 
topological correlator barely decays and it behaves in the 
same way as the staggered spin correlation function of the 
5=1/2 Heisenberg model on the square lattice. At the 
same time, the direct staggered spin correlator does not 
show any sign of order. As we already discussed, at these 
high stripe densities the 'string gas' induced stripe order 
becomes very weak and the disorderly behavior of the 
direct spin correlator is entirely due to the quantum dis- 
order in the stripe sector. Hence, Fig. ^ demonstrates 
that despite the presence of a 'hidden' spin order which is 
as strong as the Neel order found in the pure spin prob- 
lem, this can be obscured completely from the view of the 
experimentalist which can only measure the direct spin 
correlator which is affected by the fluctuating anti-phase 
boundaries. This is a lesson to keep in mind when con- 
fronted with claims that dynamical stripes do not exist 
because they do not appear in measured dynamical spin 
susceptibilities. 

VIII. DESTROYING SUBLATTICE PARITY 
ORDER. 

In the previous sections we have described a theory of 
stripe order as it emerges from a microscopy dominated 
by quantum fluctuations. In fact, when we started the re- 
search described in the above, we were after the physics of 
stripe quantum disordered states. The truely interesting 
problem is of course the phenomenon called 'dynamical 
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stripes', referring to experimental anomalies found in the 
fully developed superconductors. If it has anything to do 
with stripes, it has to be that these stripes phases are 
quantum disordered in an essential way. Hence, when 
we discovered the spontaneous directedness of the Eskes 
strings, followed by the string-gas order-out-of-disorder 
mechanism we were at first disappointed. It took a while 
to just appreciate the above on its own merit. 

We expect that any reader who is aware of the rel- 
atively developed theoretical understanding of the mi- 
croscopy of the stripes would have noticed that the 
'stripyness' as imposed by the perfect sublattice parity 
order is too literal. Starting with the true holes and 
spins as described hy t — J type models one ends up with 
a picture which is much less orderly in this regard (Za- 
anen, 1998; White and Scalapino, 1998; Morais-Smith 
et ai, 1998; Chernyshev et ai, 2000; Tchernyshyov and 
Pryadko, 2000; Martin et ai, 2000). For rcahstic values 
of parameters the holes are much more loosely bound to 
the stripes and one might even wonder if at the densi- 
ties of interest to the superconductors one can uniquely 
assign a particular real hole to a particular stripe. We 
are well aware of this and the above should be considered 
as no more than a fixed-point theory, a strong coupling 
limit with regard to the stripyness which will have a fi- 
nite basin of attraction. The physics described in the 
above will be robust against some degree of local viola- 
tion of the order, and in this restricted sense it might 
tell something about the origin of stripe order. At the 
same time, this robustness has its limitations and when 
the fiuctuations increase at some point a phase transition 
has to follow where sublattice parity order is destroyed. 
At the same time, it could well be that the local sublat- 
tice parity fluctuations disorder charge and/or spin well 
before the transition occurs where sublattice parity order 
vanishes. Therefore, a variety of distinct, partially disor- 
dered phases can exist in between the fully ordered stripe 
phase and the fully disordered state and these might have 
something to do with high Tc superconductivity (Zaanen 
and Nussinov, 2000). 

What are the fluctuations violating sublattice parity 
order? As we already stated repeatedly, perfect order of 
this kind means that stripes form trajectories of nearest- 
and next-nearest neighbor links on the lattice and to vi- 
olate this order one has to violate this connectedness re- 
quirement. The disorder excitation is therefore simple 
and unique (Zaanen and Nussinov, 2000): a stripe com- 
ing to an end, an object which we called the 'stripe dis- 
location' (Fig. |l|). 

This entity has clearly a topological status: a single 
stripe dislocations destroys sublattice parity everywhere. 
The sublattice parity in the region 'below' the stripe dis- 
location cannot be matched consistently with the parity 
'above' the dislocation (Fig. |l^ and thereby it destroys 
the notion of a definite sublattice parity. At the same 
time, it is a classic dislocation with regard to the charge 



order. It is like the half row of atoms of metallurgy which 
is well known to be the topological excitation destroy- 
ing the translational symmetry breaking associated with 
crystalline order. Finally, it represents clearly a disorder 
event with regard to Neel order - it destroys the bipar- 
titeness of the squeezed lattice and thereby it represents 
an essential spin frustration which will destroy the Neel 
order everywhere as well. However, it is not quite a gen- 
uine topological excitation of the spin system, as we will 
discuss later. 

Despite its simple appearance, this stripe dislocation 
has the remarkable meaning that it is the omnipotent 
disorder excitation belonging to stripe order. It destroys 
at the same time anti-phase boundarieness, charge- and 
spin order. Its significance becomes particularly obvious 
in combination with the general principle of duality. 

Duality is a mathematical principle with a general ap- 
plicability in continuum field theory, stating that a deep 
relation exists between states of matter separated by a 
phase transition. In the condensed matter context its 
applicability is limited to situations where one can stay 
away from the lattice (UV) cut-off, which means in prac- 
tice that it has little to say about strong first order tran- 
sitions. This relation is as follows. Given a long range 
order, excitations can be uniquely defined using the ma- 
chinery of topology which destroy this order globally. 
The field configurations of the continuum theory can be 
rigorously subdivided in smooth configurations and sin- 
gular configurations corresponding with the topological 
excitations. The smooth configurations cannot destroy 
the order because they are themselves part of the order 
and the disorder is carried entirely by the topological 
excitations. Hence, at the order-disorder transition the 
topological excitations proliferate in the vacuum. Be- 
cause these topological excitations are interacting enti- 
ties occurring at a finite density in the disordered states 
they in turn define an interacting system with a tendency 
to break symmetry spontaneously. The order parameter 
theory of this disorder matter is than equal to the field 
theory describing the long-wavelength physics of the dis- 
ordered state. We are intrigued by the following question: 
could it be that static stripes and the high Tc superconduc- 
tors are related by duality ? If this would be the case the 
physics of the superconductors should be related to the 
physics of stripe dislocation matter, because the stripe 
dislocations are the elementary disorder excitations as- 
sociated with stripe order. Different from the theory of 
stripe order, we do not understand the nature of stripe 
disorder theory at all. 

The fundamental problem is that at least for spin 
5 = 1/2 the stripe dislocations are essential frustrations 
(Zaanen and Nussinov, 2000) restoring a Marshall sign 
structure in the vacuum, and these signs represent a dif- 
ficult but not necessarily intractable problem (Weng et 
ai, 1997). Without these signs we have some understand- 
ing of the disorder theory and let us sketch some of the 
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FIG. 14. (a) The fundamental stripe dislocation which not only destroys the charge order, but also the sublattice parity 
order, while it gives rise to an essential frustration in the spin system. The bold arrow corresponds with the Burger's vector 
while the open arrows indicate the half-vortex like texture it causes in a semi-classical antiferromagnet. In the inset its action 
in squeezed space is indicated, (b) Initially a dislocation and an anti-dislocation (with regard to sublattice parity) will bind due 
to their logarithmic spin-mediated interaction. This sublattice parity 'neutral' dislocation will still destroy the charge order 
although it is no longer an essential frustration in the spin system. 



essentials. Let us assume that the charge system is made 
out of bosons while the spin system can be represented 
with the 0(3) quantum non linear sigma model, which is 
also a bosonic field theory. Except that we now allow for 
stripe dislocations, and the neglect of the Marshall signs 
of the spins, everything else is like the situation described 
in Section VII. What can happen? The answer is that 
the fully ordered stripe phase melts initially in a charge 
quantum nematic which is at the same time a spin ne- 
matic. Subsequently, this state can either undergo a first 
order transition into an isotropic quantum spin liquid, 
or in a incompressible quantum spin liquid characterized 
by a topological order. The effective field theory govern- 
ing this spin dynamics is the Z2 gauged 0(3) quantum 
non-linear sigma model. 

The qualitative idea is straightforward. Consider the 
stripe dislocation, insisting that the spin system is semi- 
classical. It is directly clear that the 0(3) order parame- 
ter will fold around the dislocation forming a half-vortex 
like texture, see Fig. |l^. Interestingly, this is not a topo- 
logical excitation in the 0(3) spin system. The topolog- 
ical excitation associated with the internal 0(3) symme- 
try in 2-1-1 D is a skyrmion (Fradkin, 1991), and this is a 
texture living on the time slice which involves rotations 
in order parameter space in two orthogonal directions. 
The half vortex associated with the stripe dislocation ro- 
tates only in one plane and it is easily seen that it carries 
in addition a zero-mode associated with the rotation in 
the orthogonal plane. Despite these intricacies, it is still 
true that on the time slice these half-vortices interact 
via long range, logarithmic interactions mediated by the 
spin system. In this sense they are like 0{2) vortices 
and these logarithmic interactions will cause dislocations 



and antidislocations to bind in pairs initially, which are 
globally equivalent to the 'neutral' dislocation indicated 
m Fig. |l|. This bound pair of dislocations is not affect- 
ing the sublattice parity globally. At the same time it is 
a dislocation with regard to the charge system carrying 
a Burger's vector which is twice that of the elementary 
dislocation. 

Referring back to section V, the order-out-of-disorder 
argument was based on the observation that dislocations 
of the type Fig. |lj cannot proliferate. However, this 
could only be proved under assumption that stripes are 
uninterrupted elastic lines. Allowing for a small but 
finite break up probability invalidates this assumption 
and thereby the argument. End points of stripes fluc- 
tuate much more than intact stripes and therefore they 
will proliferate always when the kinetic energy is large 
enough. However, because of the argument presented in 
the previous paragraph, sublattice parity order is still 
protected because of the logarithmic interactions medi- 
ated by the spin system. 

Since these neutral dislocations restore translational 
invariance, charge order is destroyed and the stripes form 
a quantum fluid. In fact, this is a nematic quantum liq- 
uid crystal, of the kind introduced by Kivelson, Fradkin 
and Emery, 1998. However, this is a subject of its own 
which is not essential to the remainder of the present 
argument. What matters is that the charge disorder im- 
plies a gap and the next question is on the nature of the 
spin dynamics at energies less than this mass-gap. 

Since the neutral dislocations do not destroy sublat- 
tice parity globally, the 'hidden' spin order of section VII 
can survive in principle. These dislocations cause local 
frustrations which will act to decrease the hidden Neel 
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order but this does not necessarily imply that the order 
completely disappears at the moment that these disloca- 
tions proliferate. Hence, in principle a state exists char- 
acterized by the hidden Neel order as measured by the 
topological correlator Eq. while the charge sector is 
quantum disordered. This state is a spin nematic. 

Spin nematics were first introduced by Andreev and 
Grishchuk, 1984, on basis of general symmetry consid- 
erations. However, to the best of our knowledge these 
spin nematics have never been identified in experiment, 
and we are claiming that this type of spin order arises 
in a most natural way in this stripe context. Consider a 
snapshot of a time-slice in Euclidean space time, of the 
kind as shown in Fig. ^. Point the finger at a particular 
site in a magnetic domain. Because the dislocations are 
proliferated, stripes are delocalized and it has to be that 
in going along the time direction at some point a stripe 
will pass this particular site and after this passage the A 
sublattice has changed in the B sublattice and vice versa. 
Therefore the Neel order parameter will point in exactly 
the opposite direction. There is still a sense of broken 
spin rotational symmetry. Call the initial direction the 
north pole. After the stripe has passed the order param- 
eter will point to the south pole. Hence, north pole and 
south pole are identified but the location of the north 
pole can still be freely chosen on the half-sphere. This is 
the director order parameter which is usually associated 
with nematic order, and therefore this state should be 
called a spin nematic. 

Traditionally (de Gennes and Prost, 1993), the ef- 
fective theory of (spin) nematics is written in terms of 
a tensor order parameter, for a three component spin 
(qap) = {nanp - l/3Sa.p) ( a,/3 = x,y,z), which is 
clearly invariant under n ^ —n (identification of the 
poles). However, it was only quite recently realized by 
Lammert, Rokshar and Toner, 1995, that the complete 
effective theory should explicitly incorporate the Ising 
gauge invariance associated with the director, which is 
automatic when the theory is written in terms of the 
redundant vector degrees of freedom n. It is an Ising 
gauge invariance because the vector is defined modulus 
its sign. For the spin-nematic, this gauge theory is as fol- 
lows in Lagrangian formulation. Define a 3d cubic lattice 
(2 space and 1 time direction) and define on every site 
a 0(3) vector fti. If these were coupled by normal ex- 
change interactions this would just correspond with the 
0(3) QNLS. However, define now the Ising variables t^j 
living on the bond between sites taking the values 
±1. The effective action describing the spin-director or- 
der parameter theory is, 

= - J ^ • fij + K^Uarf^ (28) 

<ij> □ 

where the last term is the plaquette action defining Ising 
gauge theory: take the product of the values of the r'^'s 



living on the bonds of one particular plaquette, to sum 
subsequently over all plaquette's. It is easily checked 
that this theory is invariant under the simultaneous local 
transformations Gi = J2j{i)'^ij with t^\±) = |=f), i-e. 
reversing all signs of the bonds emerging from site i, and 
the north-south identification n ^ —n. Hence, this is the 
famous Z2 lattice gauge theory, minimally coupled to an 
0(3) matter field. 

It is well known (Kogut, 1979) that the Z2 pure gauge 
theory has a confining and a deconfining phase. A rep- 
resentative gauge-fixed configuration for the deconfining 
phase is simply the one where all r's take the value -1-1. 
Adding the matter fields an keeping this configuration 
of gauge fields, an ordered state exist with all vectors 
pointing in the same direction. By repeated gauge trans- 
formations starting from this gauge-fixed configuration, 
one recovers the spin nematic order as we just discussed. 
With a little thought one can convince oneself that the 
abstract Z2 gauge degree of freedom of Eq. (|28| ) acquires 
in this stripe context a material meaning. In the ordered 
stripe phase, the sublattice parity has acts like a global 
symmetry. Having decided that a particular domain has 
sublattice parity -1-1 one can unambiguously establish the 
sublattice parity of any other domain. However, when 
charge is disordered this is no longer possible. Instead, 
at energies less than the charge gap, the sublattice par- 
ity survives as a Z2 degree of freedom governed by local 
symmetry! 

This is not all because two other phases are possible. 
The Z2 sector can stay deconfining while the 0(3) sector 
disorders (i.e., keep K large and decrease J). This cor- 
responds in the stripe language with a disordering of the 
hidden spin order, driven for instance by the frustrations 
associated with local stripe break-ups, while the sublat- 
tice parity order stays intact. Finally, the Z2 sector can 
confine. It is well understood that this confining phase 
corresponds with the ordered phase of a dual global Ising 
model where the dual variables correspond with the topo- 
logical excitations of the Ising gauge theory, the 'fluxons' 
or 'visons'. These can be visualized in the same gauge 
fix as for the deconfining phase as lines cutting through 
the bonds starting at some point and disappearing to 
infinity where all r's have reversed their sign. This is 
like a half-infinite line of antiferromagnetic bonds living 
in a ferromagnet and one directly recognizes that this is 
precise representation of the frustration caused by the 
stripe dislocations in the spin system. Hence, the con- 
fining phase of the theory Eq. ( p^ ) describes the stripe 
quantum liquid characterized by free dislocations. The 
material meaning of confinement is that at scales large 
compared to the characteristic dislocation distance the 
notion of sublattice parity no longer exists, even not lo- 
cally. 

Let us now return to the question we posed in the be- 
ginning: can the suspected hidden order associated with 
the high Tc phenomenon be related to stripes? Obvi- 
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ously, the implicit suggestion throughout this paper has 
been that sublattice parity order is a candidate which 
should be taken seriously. To be more explicit, could it 
be that the quantum-criticality of the optimal Tc su- 
perconductors is associated with the zero-temperature 
transition where the stripe dislocations unbind? Let us 
present some argiiments favoring this possibility: 

(i) We have hopefully convinced the reader that sublat- 
tice parity order is a genuine part of the vacuum struc- 
ture of the stripe phase. At the same time, it has to 
be that the very notion of sublattice parity is destroyed 
at sufHciently large doping. Although convincing experi- 
mental evidence is still missing, one would expect a more 
conventional Fermi-liquid/BCS physics if the hole den- 
sity becomes sufficiently large and sublattice parity does 
not exist as a degree of freedom in conventional fermiol- 
ogy. In the above we have spelled out the unique way in 
which sublattice parity gets destroyed. Sublattice par- 
ity can persist as a degree of freedom governed by local 
symmetry even in a state which is spin- and charge wise 
quantum disordered. This is of course the deconfining 
state of the Ising gauge theory. The meaning of the con- 
fining state is that at low energies the whole notion of 
sublattice parity has disappeared from the long wave- 
length theory. Hence, the BCS superconductor is at the 
same time a sublattice parity confining phase. Since we 
seem to know the low hole density (stripes) and high hole 
density (BCS) limits, it has to be that there is phase tran- 
sition in between where sublattice parity gets confined. 
A priori, it cannot be excluded that this phase transition 
is of the first kind where necessarily the disappearance of 
sublattice parity goes hand in hand with other symme- 
try changes like charge- and/or spin disordering (e.g., the 
first order transition in the spin nematic where both n 
gets disordered and Z2 confines). However, this appears 
as unlikely, given the abundant evidence for quantum- 
criticality associated with the stripes. 

(ii) Above all, sublattice parity is a very hidden degree 
of freedom. There is no existing experiment which can 
directly measure if sublattice parity order exists or not, 
especially so in the charge- and spin disordered phases 
where it only exists as a local degree of freedom. It does 
have indirect consequences which are accessible to ex- 
periment at least in principle. For instance, it can be 
easily seen that the elementary excitations of the quan- 
tum spin-nematic are associated with S = 2 instead of 
the usual triplets (Zaanen, unpublished) and these can- 
not be directly measured with neutrons (see also Andreev 
and Grishchuk, 1984). It could well be that massive ex- 
citations of this kind exist because they are again very 
hard to measure. 

(iii) There has been little mention of superconductivity. 
The reason is that superconductivity comes for free, at 
least under the assumption that the stripe phase itself is 
associated with pairs of electrons. Charge- wise the stripe 
phase is then a bosonic crystal and superconductivity 



cannot be avoided when the charge quantum disorders. 
Just as in the case of superconductivity and antiferro- 
magnetism (Zaanen, 1999; van Duin and Zaanen, 2000) 
the mode-couplings between the superconducting phase 
mode and all other modes are supposed to be irrelevant 
operators and therefore superconductivity can be simply 
superimposed on the physics of the spin- and sublattice 
parity sectors, at least at zero temperature, 
(iv) There is one very serious problem: there is no obvi- 
ous place for the nodal fermions in this framework. This 
is partly by construction. In the above we have assumed 
that everything is bosonic and in a bosonic universe there 
is no room for S* = 1/2 excitations. Nodal fermions have 
surely to do with this spin quantum number. We already 
announced that we made crucial assumption by neglect- 
ing the Marshall sign's. As is well known, a S = 1/2 
Heisenberg spin problem defined on a bipartite, nearest- 
neighbor bond lattice has a bosonic ground state in the 
sense that the ground state wave function is nodeless. 
The reader familiar with this construction will immedi- 
ately infer that in the presence of the stripe dislocations 
this is no longer true. Hence, at the moment stripes are 
no longer perfectly connected the spin system acquires a 
non-trivial sign structure in vacuum and one no longer 
knows what to expect (Weng et al., 1997). Is this the 
missing ingredient, linking the above to high Tc super- 
conductivity? We have no clue, but we hope to have con- 
vinced the reader that it is very unreasonable to believe 
that the relationship between stripes and superconduc- 
tivity is understood. 

Note added in press: after completion of this 
manuscript we learned that the prediction of Zaanen, 
2000, of an induced modulus string gas which has a 
stretched-exponential dependence on the density has 
been confirmed by numerical simulations (Yoshihiro 
Nishiyama, preprint, Okayama University) . 
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